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ABSTRACT 

We present new radial velocity results for 176 stars in the Fornax dwarf spheroidal galaxy, of which at least 
156 are probable Fornax members. We combine with previously published data to obtain a radial velocity 
sample with 206 stars, of which at least 176 are probable Fornax members. We detect the hint of rotation 
about an axis near Fornax's morphological minor axis, although the significance of the rotation signal in the 
galactic rest frame is sensitive to the adopted value of Fornax's proper motion. Regardless, the observed stellar 
kinematics are dominated by random motions, and we do not find kinematic evidence of tidal disruption. The 
projected velocity dispersion profile of the binned data set remains flat over the sampled region, which reaches 
a maximum angular radius of 65'. Single-component King models in which mass follows light fail to reproduce 
the observed flatness of the velocity dispersion profile. Two-component (luminous plus dark matter) models 
can reproduce the data, provided that the dark component extends sufficiently beyond the luminous component 
and the central dark matter density is of the same order as the central luminous density. These requirements 
suggest a more massive, darker Fornax than standard core-fitting analyses have previously concluded, with 
M jLy over the sampled region reaching 10 to 40 times the MjLy of the luminous component. We also apply 
a non-parametric mass estimation technique, introduced in a companion paper Although it is designed to 
operate on data sets containing velocities for >1000 stars, the estimation yields preliminary results suggesting 
MjLy ^ 15 inside r <1.5 kpc. 

Subject headings: galaxies: dwarf — galaxies: individual (Fornax) — galaxies: kinematics and dynamics — 
(galaxies:) Local Group — methods: statistical — techniques: radial velocities 



1. INTRODUCTION 

The Milky Way's dwarf spheroidal (dSph) satellite galax- 
ies have stellar masses similar to those of globular clusters 
(Miummous ^ IQ^~''Mq), yet they are much more spatially ex- 
tended (R - 0.5-3 kpc for dSphs; R - 0.01-0.05 kpc for 
globular clusters). These characteristics give dSphs the low- 
est luminosity densities of any known galaxies. The discov- 
ery that their internal velocity dispersions all exceed 7 km 
s~' (Aaronson 1983; Mateo 1998 and references therein) has 
given rise to competing interpretations and speculations con- 
cerning their origin and cosmological significance. 

If dSphs are assumed to be approximately virialized sys- 
tems, their large velocity dispersions indicate the presence 
of copious amounts of dark matter Estimates of mass-to- 
light ratios (M/L) based on the equilibrium assumption have 
yielded M/L ^ 5-500 in solar units for various Milky Way 
dSphs (Mateo 1998 and references therein; Kleyna et al. 
2001; Odenkirchen et al. 2001; Kleyna et al. 2005). The 
dSphs are then the smallest, nearest systems believed to re- 
side within dark matter halos, and so provide a convenient 
and fundamental testing ground for cold dark matter models 
of structure formation. 

Alternatively, large measured velocity dispersions have 
been cited as possible evidence that the dSphs are presently 
undergoing tidal disruption as they orbit within the MW po- 
tential (Kuhn 1993; Kroupa 1997; Klessen & Kroupa 1998; 
Klessen & Zhao 2002; Fleck & Kuhn 2003). According to 
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this interpretation, dSphs may be in a state far from dynam- 
ical equilibrium, and masses derived under that assumption 
may be inflated. If the observed stellar velocity dispersions 
can be attributed to streaming tidal debris projected along the 
line of sight, the need to invoke dark matter for explaining 
dSph kinematics subsides, perhaps entirely (Fleck & Kuhn 
2003; Kroupa 1997). 

Large samples of radial velocities measured for dSph stars 
may be capable of distinguishing between various equilibrium 
and tidal models by examining the velocity trends across the 
face of the system. If a dSph is close to dynamical equilib- 
rium, its stellar motions provide an estimate of the underlying 
mass distribution. Tidal disruption is expected to be accompa- 
nied by a radial velocity gradient, giving rise to apparent rota- 
tion with a characteristic orientation. We present in this paper 
new radial velocity results for 176 stars along the line of sight 
to Fornax. After combining with previously published results, 
we test for rotation and then measure the radial velocity dis- 
persion profile extending from the Fornax center to the nomi- 
nal edge of the luminous component. We consider the results 
in the contexts of equilibrium and tidal disruption models. We 
also estimate the Fornax mass non-parametrically, applying a 
technique formally introduced in a companion paper by Wang 
et al. (2005; hereafter. Paper I). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Photometry, Astrometry, and Target Identification 

In order to identify spectroscopic target stars, we first ob- 
tained photometric data from 31 fields covering a 110' x 90' 
region of sky centered on the Fornax dSph (q!2ooo=2:39:52, 
^200o=-34:28:09). These observations took place the nights 
of 1993 November 30, December 1, and December 10, dur- 
ing photometric conditions. The data consist of 600-700s ex- 
posures in both V and I filters using flie 2048 x 2048 TEK 
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3 CCD detector at the Las Campanas Observatory 40-inch 
telescope (field size=24'x24', scale=0.7 arcsec pixel"'). The 
images were processed using twilight flatfield exposures and 
multiple bias frames. We used the two-dimensional stellar 
photometry program DoPHOT (Schechter et al. 1993) for the 
reductions and placed the resulting instrumental magnitudes 
on the Kron-Cousins scale (Bessell 1976) using 56 Landoldt 
photometric standard stars observed during the same nights. 
From the formal error values returned by DoPHOT, and mul- 
tiple measurements of stars in overlapping field regions, we 
estimate our photometric accuracy to be ±3%. 

The resulting color-magnitude diagram (CMD) is shown in 
Figure 1 , which also shows the region in the CMD from which 
we selected stars for spectroscopic observation. The bound- 
aries of this region were chosen to enclose points representing 
likely Fornax members, as well as those stars bright enough 
to maintain reasonable integration times during spectroscopy. 
The chosen region roughly extends along the brightest ~ 1 .5 

1- band magnitudes of the Fornax red gianl branch (RGB) and 
includes more than 4000 Fornax candidate members. We 
conservatively chose from among the brighter targets in this 
selection region for spectroscopic foUow-up. Based on our 
spectroscopy results, it is feasible to derive accurate radial ve- 
locities for stars at even the faint edge of this selection region. 

Since one of our spectroscopic observing runs took place 
prior to 1993, some targets were selected without the bene- 
fit of this photometry. For target identification leading up to 
the 1992 November-December spectroscopy run, we relied on 
photometry obtained in 1990 November over a smaller region 
of sky. Reduction of these photometric data is described in 
Mateo et al. (1991; hereafter, M91). Many of these targets fall 
outside the selection region shown in Figure lb. This earlier 
effort was devoted to observing primarily the brightest candi- 
date members, although several faint RGB stars were selected 
in order to probe the limits of the instrumentation. 

An additional factor entering our target selection was a 
star's sky position relative to the center of Fornax. To aid this 
selection we converted the (x,y) CCD position returned by 
the DoPHOT centroid algorithm into equatorial coordinates 
using the IRAF routines TFINDER and CCTRANS, and tied 
our astrometry to the USNO-IB system using up to several 
hundred USNO stars per CCD frame. From measurements of 
stars in overlapping fields, we estimate the 2cr astrometric ac- 
curacy to be better than 0.2". In deciding on eventual spectro- 
scopic targets, a selection routine closely following the stellar 
density distribution is inadequate. The outer, sparsely popu- 
lated regions are of disproportionately high kinematic inter- 
est. Nevertheless, we wished to obtain a large sample with a 
high fraction of Fornax members. In the end we chose at least 

2- 3 candidate members in all the outer Fornax CCD fields, 
and limited our selection in the iimer fields to ~ 5-8 stars per 
field. Figure 2 maps the locations of the stars falling within 
our CMD selection region, and identifies which of those stars 
we ultimately observed spectroscopically. 

2.2. Spectroscopy 

We obtained spectra of the selected red giants over four 
observing runs occurring 1992 November 28 - December 7, 
1993 December 12-20, 1994 October 22-28, and 2002 De- 
cember 12-15. The 1992-1994 spectra were taken at the 
Las Campanas Observatory 2.5m du Pont Telescope equipped 
with an echelle spectrograph and 2D-Frutti photon counting 
detector (Shectman 1984). The 2002 spectra were acquired at 
the Magellan 6.5m Clay Telescope with the Magellan Inamori 
Kyocera Echelle (MIKE) spectrograph (Bernstein et al. 2003) 



using a 1.0" x 5.0" sHt and MIKE's red-side CCD detector 
set to obtain spectra between 4400-8000 A. As they came 
from different telescopes/instruments, the 1992-1994 spectra 
and 2002 spectra were reduced independently. We followed 
the same general procedure in both cases. In what follows, 
we describe the observation and reduction procedure specific 
to the 2002 spectra. Where details differ regarding the 1992- 
1994 spectra, we comment parenthetically. 

Object exposure times ranged from 360-720 seconds (900- 
2400s for the 1992-1994 runs), with most at 600s (1800s). In 
addition, we took Is (360s) exposures of a Th-Ar compari- 
son arc lamp before and after object exposures. Other expo- 
sures included quartz-illuminated dome flats and spectra of 
bright radial velocity standard stars used for aperture identi- 
fication and velocity calibration. We used the IRAF data re- 
duction software package to reduce the raw spectra to helio- 
centric radial velocities. After overscan and bias subtraction, 
we produced a master flatfield frame for each night by averag- 
ing quartz-illuminated dome exposures. Flatfield frames were 
then normalized using the IRAF task APFLATTEN, which 
models and removes both the spatial profile and spectral shape 
of the illumination pattern, leaving only the sensitivity varia- 
tions. We corrected our object spectra for these sensitivity 
variations by dividing each by the appropriate normalized flat- 
field frame. We then ran the IRAF task APALL (with FOR- 
MAT keyword set to "strip") on the spectra of the bright stan- 
dard stars to obtain two-dimensional traces of the echelle or- 
ders on the detector and rectify the spectra. We found the 
order location on the detector to remain quite stable over the 
course of the run (MIKE is fixed with respect to gravity). We 
then ran APALL on all the target object spectra to obtain a rec- 
tified spectrum from each order, referencing the trace pattern 
identified for the most recently observed bright star. Thorium- 
argon comparison spectra were rectified in exactly the same 
manner as the individual stellar spectra they would eventually 
correct and calibrate. 

For the 2002 data we had also to remove cosmic rays (the 
1992-1994 data obtained via photon-counting device did not 
suffer from cosmic rays). As a first pass, we ran the IRAF 
task COSMICRAYS, specifying conservative thresholds so 
as to remove only the most conspicuous events. The major- 
ity of cosmic rays were then removed by the task CONTIN- 
UUM, using a lOth-order cubic spline to replace any pixel 
value above an upper sigma threshold determined by eye to 
optimize accurate cosmic ray identification. 

We then employed the tasks IDENTIFY and REIDENTIFY 
to convert the Th-Ar spectra from pixels to wavelength space. 
A typical arc lamp spectrum would have, for the ten orders 
(four orders) we eventually used for velocity measurement, 
140 lines (320 lines) reidentified with an rms scatter of 0.09 
A(0.03 A) for a fourth-order polynomial fit to the wavelength 
solution. These wavelength calibrations were then applied to 
the object spectra using the DISPCOR task, which converts 
units in the dispersion direction from pixels to Angstroms us- 
ing the weighted solutions to the two nearest comparison arc 
spectra as references. 

Up to this point, the spectra remained two-dimensional, and 
we had treated each row in the spatial direction separately. 
This is necessary because the spatial and spectral axes are not 
orthogonal within each order (i.e., the spectral lines are tilted 
by ~ 20° with respect to the spatial direction on the detector). 
By wavelength-calibrating each row in the spatial direction 
explicitly, we eliminate this problem and retain the full spec- 
tral resolution of the instrument. 

We then converted the wavelength-calibrated two- 
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Fig. 1 . — Fornax red giant branch, a) includes all stars measured photometrically in a 1 10' X 90' region of sky centered on the Fornax dSph. b) shows only those stars observed 
spectroscopically. and illustrates the boundaries of our color-magnitude target selection region. Filled circles are probable Fornax members, based on velocity criteria described in 
section 3.2. Open circles are probable foreground contaminants. Open triangles represent stars with marginal membership status. Points located outside the CMD selection region 
represent stars observed for this study before the photometry comprising this CMD was available, and so were chosen based on the photometry described in Mateo et al. (1991). 
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Fig. 2. — Maps of (a) all stars meeting the selection criteria discussed in Section 3.2; overplotted are the boundaries of the 31 photometric fields observed, (b) maps stars for which 
we measured radial velocities. Filled circles represent stars later determined to be probable Fornax members. Open circles represent stars rejected as probable foreground contaminants 
on the basis of their radial velocities. Open triangles represent stars having velocities marginally consistent with Fornax membership (see section 3.2). The inner and outer ellipses 
are the King core and tidal radii, respectively, which have pubUshed semi-major axis values rc„„ = 13.7' ± 1.2' and r,ijg = 71.1' ± 4.0', with ellipticity e = 0.30 ± 0.01 (Irwin & 
Hatzidimitriou 1995). The standard coordinate system is centered on the Fornax dSph such that ({,?;)=(0,0) corresponds to a20oo=2:39:52, 52ooo=-34:28:09. North is toward the top of 
the figure, east is to the left. 



dimensional spectra to one-dimensional spectra (for the 
1992-1994 spectra, this was accomplished by the ECDIS- 
PCOR subroutine). For the 2002 spectra, first we used 
SCOMBINE to sum the five spatial rows at the center of each 
spectral order, which we had determined to carry the stellar 
signal. Separately, we averaged two rows located sufficiently 
far from that aperture center so as to identify primarily the 
sky spectrum. Finally, we used SCOMBINE to subtract the 
normalized sky spectrum from the summed stellar spectrum 
of the same aperture. A second pass with CONTINUUM then 
fitted and subtracted the continuum signal, and IMREPLACE 



then replaced with zero any pixel with absolute value greater 
than fifty. Our velocity measurements are limited to the 

wavelength range - 4900-5600 A(- 5120-5460 A), to im- 
prove precision of the results. The wavelength range chosen 
for the 1992-1994 data corresponds to the four echelle orders 
determined to yield the most precise velocity measurements 
from these spectra. At wavelengths redder than this range, 
the metal-poor spectra provide few absorption lines, while at 
wavelengths blueward of 4900 A, there is little continuum 
flux from the red giant targets. The strongest absorption 
lines within the selected region belong to the magnesium-B 
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triplet having a rest wavelength near 5170 A, while the many 
weaker absorption lines contribute usefully in aggregate to 
the cross correlation function. 

We calculated heliocentric radial velocities using the cross- 
correlation package FXCOR. We cross-correlated the ex- 
tracted spectrum of each target star against a high signal-to- 
noise ratio template consisting of the sum of 27 (75) spectra 
of bright radial- velocity standard stars. The spectrum for each 
standard had first been shifted to a common heliocentric red- 
shift equal to that of the star HD6655 (vte/,o=19.5 ± 0.3 km 
s"^; Udry et al. 1999). Prior to cross-correlation, we filtered 
the 10970 (8192)-pixel spectra with a ramp function, cutting 
on at 175 (100) wavenumbers, increasing linearly to full value 
at 200 (170) and then decreasing linearly from 2200 (700) to 
a cutoff at 2500 (1000). A Gaussian fit then located the cen- 
ter of the cross-correlation peak, thereby specifying the radial 
velocity difference between the object and template. 

2.3. Measurement Uncertainties 

As a check on the reliability of the extractions and cross- 
correlations, we compare independent velocity results ob- 
tained for the bright standard stars observed multiple times. 
Let Nbright be the number of standard stars observed, and let 
Nb be the number of independent observations of standard 
star b. Letting Nb be the total number of individual mea- 
surements accumulated for standard stars {Nb = Y^h=i'" ^b), 
we then define the cumulative sample variance over all in- 
dependent measurements of such stars to be crlright - (^b - 
Nbrightr'Y:tr'T,'J=M.j-{v)bf. For our sample we find 

(^bright = 0.89 km s"' iabright = 0-72 km s"') for Nb=24 
(Nb=107). This indicates a satisfactory level of internal con- 
sistency for our purposes. 

To calculate the internal measurement uncertainty, aj, as- 
sociated with each independent velocity measurement, vj, we 
assume that multiple measurements of a given star having 
true velocity Virue follow a Gaussian distribution with mean 
Vtrue and variance aj. Multiple measurements will be dis- 
tributed as Vtrue + where ej is a random variable fixed 
by measurement and following a Gaussian probabiUty distri- 
bution with mean zero and variance unity (a standard nor- 
mal distribution). For simplicity we estimate v,,.„,..,, the true 
velocity of star /, from A',- independent measurements as^ 
Vtrue,i = 12%\ ^ij- ^c make the further assumption that 
the difference v,„,e., - v,„,g., is negligible. This assumption is 
perhaps naive; however, we find for this data set that a rig- 
orous treatment, properly considering the uncertainty in Vtrue, 
gives nearly identical results. With these assumptions we then 
express the velocity measurement of star i as 

Vij = Vtrue,i + <7ij<^ij, (1) 

We model the cr,; as a function of the Tonry-Davis R value 
(Tonry & Davis 1979), which FXCOR calculates as the ratio 
of the selected cross-correlation peak height to the average 
height of the nonselected peaks. We express the relationship 
as 



This two-parameter model generalizes the original Tonry & 
Davis formalism, which assumed x= I. We find that we are 

^ We follow the convention by which the estimation of quantity q is de- 
noted g. 



better able to reproduce the empirical sample variances ob- 
tained from repeat measurements by treating x as a free pa- 
rameter. 

The base- 10 logarithm of the squared error in the i'^ mea- 
surement is then 

log[iVij-v,rue,if] = 21ogQ-2xlog(l +Rij)+logiejj). (3) 

The term log(e? ) has mean (log(e^y)) = -0.55, from Monte 
Carlo simulations. If we define 6ij = log(ey) + 0.55. Equation 
3 becomes 

log[(v,;-v,™,,,)^] = 21oga-2xlog(l +/?,;■) + 0.55. (4) 

We then estimate x and a by linear regression using the 
{Rij,Vij) data from only those stars with repeat measurements, 
and recognizing that the Sij have mean value {Sij) = 0. In- 
cluding bright standard stars, our 2002 data contain 25 (139 
for the 1992-1994 data) repeat observations of 6 (19) differ- 
ent stars. From the 1992-1994 data we obtain the estimates 
a = 6.0 km s~' and x = 0.50. Because the 2002 observations 
contain fewer repeats, and no repeat observations of low-/? 
target stars, we adopt x = 0.50 for the 2002 data and then es- 
timate a with a least-squares fit to find a = 7.6 km s~' for the 
2002 results. Using the appropriate values for the parameters 
a and x, we then calculate the uncertainty in each individual 
velocity measurement using Equation 2. 

Within the 1992-1994 data, we were able to check the sta- 
bility of our velocity zero-point both night-to-night and run- 
to-run using individual spectra obtained each night for the 
bright star CPD-35°919, located just west of the Fornax cen- 
ter. Within a given run, the night-to-night scatter in the mea- 
sured velocity of CPD-35°919 is nearly identical to its esti- 
mated internal errors. Comparing the mean velocity measured 
for this star in each run, we find no significant run-to-run scat- 
ter. 

Since the 1992-1994 and 2002 spectra were cross- 
correlated using different template spectra, we searched for 
any systematic zero-point velocity difference between the 
two independent data sets. To accomplish this, we cross- 
correlated the target spectra from the 2002 run against the ra- 
dial velocity template used for the 1992-1994 data. We found 
the resulting measured velocities to differ by less than 0.05 
km s~^ with respect to their values derived from the 1992- 
1994 velocity template. Finding an equally small discrepancy 
when cross-correlating the 1992-1994 target spectra against 
the 2002 template, we take any zero-point velocity offset be- 
tween the two data sets to be insignificant. To measure the 
overall zero-point offset of the entire combined data set, we 
considered the measured velocities of all bright radial veloc- 
ity standard stars separately, finding a mean discrepancy of 
{v observed -V published) =1-25 km s"', which wc take to be our 
zero-point error. Subtracting this value from all our velocity 
results, we then measure a mean value of -0.3 ±1.2 km s~^ 
for eight twilight sky spectra (corrected individually for diur- 
nal and annual motions). Thus we are confident that we have 
placed the velocities on a true heliocentric zero point. 

3. RESULTS 

3.1. Heliocentric Radial Velocities 

In Table 1 (tables are located after the appendix) we present 
the heliocentric radial velocities and formal uncertainties de- 
rived from each individual observation. Entries are sorted by 
date of first observation, with any repeat measurements of 
the same star fisted directly below. Additional information 
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includes the equatorial coordinates, date and time of obser- 
vation, and the measured I,V -I photometry for each object. 
The distance R is the angular distance between the center of 
Fornax, taken from Irwin & Hatzidimitriou (1995; hereafter, 
IH95) to be a2ooo=02:39:52.3, 520oo=-34:28:09.0, and the pro- 
jection of the radial position vector on the plane centered on 
these coordinates. The position angles are defined to be 0° 
due north at the tangent point and 90° due east. The final 
column indicates whether we judge the star to be a proba- 
ble member ("Y") or nonmember ("N") of Fornax, based on 
photometric and velocity criteria (section 3.2). Several stars 
present borderline cases for membership, and we mark their 
membership status as "?". If the membership status is other 
than "Y," the superscript indicates whether this is due to the 
star's photometry ("p"), velocity ("v"), or both ("v,p"). We 
include all Fornax targets, as well as the bright star CPD- 
35°919. Table 1 Usts the radial velocity results for the ob- 
served standards. 

We observed several stars on multiple occasions. In sub- 
sequent analyses we take the heliocentric radial velocity of 
each multiply-observed star to be the average of that star's 
individual velocity measurements weighted by their respec- 
tive uncertainites. Table 1 gives the weighted mean velocity, 
value, and the probability, p(x^), for stars having multiple 
measurements. Several of the Fornax stars in our sample have 
velocities previously published by M91. Table 1 compares 
our measurements for these stars with those of M9I. For 1 1 
of the 14 stars common to both data sets, we find agreement 
to within the measurement uncertainties. Of the remaining 
three, two have velocities reported by M91 differing from our 
measurements by ^ 2.5(t, while the third differs by ^ IQa. 
See section 5 for a discussion on velocity variabiUty and its 
effects on our results. 

3.2. Fornax Membership 

We identify and exclude from our sample those stars that 
are likely to be foreground contaminants. Having passed posi- 
tional and photometric criteria, these interlopers are best iden- 
tified as outliers in the observed velocity distribution. The he- 
liocentric radial velocity of Fornax is ~53 km s~\ so the ve- 
locity distribution of its stars overlaps that of foreground stars 
near v ~ km s~' . This is apparent in Figure 3a, which depicts 
the distribution of the radial velocities listed in Table 1 . As 
the derived mass of a pressure-supported system scales as the 
square of velocity dispersion, it is imperative that we obtain 
a sample with minimal contamination from non-members. In 
order to accomplish this objectively, we adopt the robust bi- 
weight estimator (Beers et al. 1990) which determines a char- 
acteristic distribution width, atw, equal to the standard devi- 
ation in the special case of a normal distribution. Since 99% 
of the members in a normally distributed sample are located 
within 2.58iT of the mean, we select as a membership criterion 
I v, - (v) I < 2.5Eai,H,, and iteratively remove those stars failing 
to satisfy this condition. This rejection process converges af- 
ter four iterations, identifying 20 stars as probable foreground. 
We are left with a sample of 156 new stars we consider to be 
members of the Fornax dSph. Their radial velocity distribu- 
tion is shown in Figure 3b. 

On examination of Figure 3a, which specifies the iteration 
that removes each of the rejected stars, one may reasonably 
wonder whether iterations 4, 3, and possibly 2 of the rejec- 
tion algorithm remove what are actually Fornax member stars. 
The eye is tempted to include the three stars with radial veloc- 
ities in the range 82.5-92.5 km s~' in the wing of the distribu- 
tion centered on Fornax's systemic velocity. Reinstating these 



stars as probable Fornax members and then forcing symmetry 
on the overall Fornax distribution would argue for the addi- 
tional reinstatement of the rejected stars falling in the 12.5- 
22.5 km s~^ range. Where it is practical, we examine the ef- 
fects on our results of retaining stars rejected in the second, 
third and fourth iterations, and note that the true member- 
ship of Fornax probably includes some, but perhaps not all, 
of these stars. Pending a larger data set, we leave the mem- 
bership status of these borderline cases an open question. 

Hereafter we combine our new velocity results with the pre- 
viously published sample of M91. The 44 stellar velocities 
measured by M91 were drawn from stars belonging to one 
of two distinct fields: one centered on the Fornax core, an- 
other located along the major axis 25 arcmin southwest of 
the center. We recalculate weighted mean velocities for any 
stars measured multiple times and/or in both data sets, using 
the quoted uncertainties for each observation. In Figure 3c 
we show the radial velocity distribution of this combined data 
set, now consisting of velocities for 206 stars. The algorithm 
described above for membership determination then rejects 
(again converging after four iterations) 28 stars as probable 
foreground, including all 20 of the stars that had been rejected 
before the addition of the M91 data. The eight additional re- 
jected stars come from M91 alone, and were rejected in that 
study as well. Figure 3d shows the velocity distribution of the 
176 stars retained as probable Fornax members. Again, we 
examine the effect of retaining those stars rejected in itera- 
tions three and four (giving a sample with N=182 members), 
as well as iterations two, three, and four (giving N=186 mem- 
bers). 

3.3. Rotation and the Proper Motion of Fornax 

The relative motion between the Sun and Fornax con- 
tributes a velocity component, Vrd(l,b), along the line of sight 
to each Fornax star. A given star's galactic rest frame (GRF) 
radial velocity, Vr,GRF, is related to its heliocentric rest frame 
(HRF) velocity^, Vr^HRF, by 

Vr,GRF = VrMRF + Vrel(J ,b). (5) 

A gradient in Vreiihb) across the face of an object as large in 
solid angle as Fornax will tend to produce a spurious gradient 
in the HRF radial velocities (see Lin & Dong 2006, in prep). 
A non-rotating object might thereby give the appearance of 
rotation to the HRF observer, and a truly rotating object may 
appear to rotate at a different speed and/or about a different 
axis. In order to test for Fornax rotation we correct for this 
perspective effect by placing our HRF radial velocity data set 
in the GRF. 

Let v» be the GRF space velocity of a given Fornax star; 
we seek to determine the component of v* along the line of 
sight from the Sun to the star. Let be the velocity of the 
Sun with respect to the LSR, and let Vf be the bulk velocity 
of Fornax with respect to the LSR. Then the projection of the 
relative motion between the Sun and Fornax along the line of 
sight to the given star is the sum of scalar products: 

Vrelihb) = -. r • (VQ-Vf). (6) 

|v,| 

We apply Equations 5 and 6 along the line of sight to every 
star in the Fornax velocity data set in order to determine each 

* In either rest frame, a star's "radial" velocity is the velocity along the 
line of sight from the Sun to the star. GRF radial velocities are defined as 
those measured by an observer at the Sim's location, but at rest with respect 
to Fornax. 
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Fig. 3. — Heliocentric radial Velocity distribution of a) all 176 Fornax candidate member stars whose velocities are presented in Table 1. Stars that were later rejected by an iterative 
membership determination algorithm are numbered according to which iteration rejected them (e.g., "l"^; 1st iteration); b) the 156 stars determined to be probable Fornax members; c) 
all 209 Fornax candidate member stars after combining our data with that of M91. Again, numbers specify which iteration removed probable nonmembers; d) the 176 probable Fornax 
members from the combined data set. In (b) and (d), a thick vertical line marks the mean velocity of members calculated using maximum UkeUhood statistics. The regions enclosed by 
dotted lines in (b) and (d) represent those stars rejected in iterations 2,3 and 4. We consider these to be borderUne members. 



Star's GRF radial velocity. For each star's GRF radial velocity 
uncertainty we adopt the corresponding measurement uncer- 
tainty originally estimated in Section 2.3. We adopt the values 
Vq = 13.7 km s-i toward (/,fo)=(26.6°,31.4°) (Dehnen & Bin- 
ney 1998). The three components of Fornax's velocity with 
respect to the LSR are computed from Fornax's heliocentric 
radial velocity and proper motion, via Equations 44 - 46 of Pi- 
atek et al. (2002). We adopt +53.3 km s~' as the heliocentric 
radial velocity of Fornax. Piatek et al. (2002) and Dinescu 
et al. (2004) provide independent measurements of the For- 
nax proper motion. Since their results agree only at the ~ 2cr 
level we consider both cases independently, giving two possi- 
ble GRF radial velocity data sets. 

We then test the two resulting GRF radial velocity data sets 
for rotation. In both cases we consider the position angle of 
every star in the intermediate N=182 member sample to coin- 
cide with that of a prospective rotation axis. For each star we 



bisect the face of Fornax with a line having that star's position 
angle, and then calculate the mean GRF radial velocity from 
the member stars on either side of the line. Figure 4 plots 
the hemispheric mean velocity difference as a function of the 
bisecting hne's position angle. Panels (b) and (c) of Figure 
4 depict the GRF rotation signal assuming the Piatek et al. 
and Dinescu et al. proper motions, respectively. For compar- 
ison, panel (a) shows the HRF apparent rotation signal, un- 
corrected for any perspective-induced velocity gradient. The 
half-amphtude of the sinusoid fit in each plot measures a char- 
actersistic rotation speed, whereas the sinusoid's phase indi- 
cates the orientation of the rotation axis. We summarize the 
results of this test in Table 5. Columns 4 and 5 of Table 5 list 
the characteristic rotation speed and orientation of the rotation 
axis, respectively. The uncertainties given for these values re- 
flect the range of values obtained using all proper motions al- 
lowed within the pubhshed (la) proper motion uncertainties. 
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Uncertainties in the solar motion and in the HRF radial ve- 
locities of the Fornax stars are not considered here. To assess 
the significance of a rotation detection, we performed Monte 
Carlo simulations in which 10'' samples of 182 stars having 
positions of those in the actual sample were drawn at ran- 
dom from a non-rotating, Gaussian velocity distribution with 
a= 12.4 km s"' . Column 6 of Table 5 gives the percentage of 
these artificial samples for which we would measure a rota- 
tion speed greater than the observed speed listed in column 4. 
A lower percentage indicates a more statistically significant 
observed rotation. 

We find that although both published proper motion mea- 
surements imply Fornax rotation about an axis at position an- 
gle ~ 115°, only the rotation detected using the Piatek et al. 
(2002) proper motion is (marginally) statistically significant. 
Fewer than 10% of Monte Carlo trials produce rotation as 
fast as the ^2.5 km s"' implied by the Piatek et al. (2002) 
proper motion. Nearly three in four trials produce the ^1.2 
km s"' rotation implied by the Dinescu et al. (2004) proper 
motion. Therefore, due to perspective effects and the existing 
uncertainty in the proper motion, we cannot state definitively 
how or even if Fornax rotates. If we simply assume that For- 
nax does not rotate, we can use the apparent rotation seen in 
the HRF radial velocity data indirectly to "measure" Fornax's 
proper motion. The uncorrected HRF data indicate ~ 2.0 km 
s"^ "rotation" about an axis at ~ 140°. A Fornax GRF proper 
motion of (/i/,/ii)=(-52,+41)(mas century"'), when applied 
to these data, would produce a GRF radial velocity data set 
showing zero rotation (panel (d) of Figure 4). 

These results and Fornax's velocity dispersion of > 10 km 
s~' (Section 3.4) indicate that, aside from a possible tidal in- 
terpretation (Section 4.3), any real rotational component is dy- 
namically insignificant. Given the proper motion ambiguity, 
we use the HRF radial velocity values of Table 1 in the ve- 
locity dispersion calculations that follow. We demonstrate in 
Section 3.4 that velocity dispersions measured in the HRF dif- 
fer negligibly from their plausible GRF values. 

3.4. Velocity Dispersion 

3.4.1. Maximum Likelihood 

We use maximum likelihood statistics to estimate the mean 
heliocentric velocity and intrinsic velocity dispersion of those 
stars we have determined to be members. Let v„ and cr, 
now be the measured radial velocity, the true radial velocity, 
and the internal measurement uncertainty, respectively, for the 
i"' of A'^ stars. Then v, = Ui + aie,, where the values {ei, e^} 
have a standard normal probability distribution. There are two 
sources of variability in v, : the random, internal measurement 
uncertainty, ct,, and the intrinsic radial velocity dispersion, de- 
noted (7p, for the stars in the sample. The latter is the physical 
quantity of interest. If we assume that the values {vi , ...,vn} 
have a Gaussian distribution centered on the mean true ve- 
locity, denoted (m), then their joint probability function is the 
product of their individual Gaussian probabilities: 



P({vi 



,vn 



})=n 



1 



1=1 



:exp 



l (v,- -(«))' 
2 (af + aj) 



(7) 



Estimates of (m) and (Tp, denoted (m) and dp, are determined 
numerically as the values that maximize the natural logarithm 
of the probability function. 



Hp): 



1 



2^ 

!=1 



Vln((jf + a2) 



1^ (v,-(m))^ N 



As the logarithm is a monotonic function, this is equivalent 
to maximizing p itself (Rice 1995). To estimate confidence 
intervals for (m) and dp we recognize that the Gaussian prob- 
ability distributions for ((m) - («)) and (dp-ap) have centers at 
zero and a joint variabihty described by a covariance matrix. 
This covariance matrix. A, has elements 



A = 



a c 
c b 



(9) 



where diagonal elements a and b are the variances of {u) and 
<jp, respectively. We determine a and b from the inverse of the 
covariance matrix, which has the property 



9'ln(p) 
'WF 



d{u)d(Tp 



dapd{u) 



(10) 



Let Za/2 denote the § quantile of the standard normal distri- 
bution. For confidence intervals containing the physical val- 
ues (m) and (Tp with 100(1 -q:)% probability, we report the 

mean velocity and velocity dispersion as (m) ±Za/2y/a and 
dpiZaftVb, respectively. Conventional 68% confidence in- 
tervals are given by (m) ± ^/a and dp ± y/b. 

We estimate the mean true velocity and intrinsic velocity 
dispersion along the line of sight for three successively less 
stringent levels of membership discrimination. In case (a) we 
consider as Fornax members only the N=176 stars surviving 
all four iterations of the velocity rejection algorithm; in case 
(b) we reinstate the six stars rejected by iterations 3 and 4; 
in case (c) we further add the four stars rejected by iteration 

2. With 68% confidence intervals about the ({u),dp) pairs 
in units of km s~', case (a) gives global values (53.3 ± 0.8, 
11.1 ±0.6); case (b) gives (53.0 ±0.9, 12.4 ±0.8); case (c) 
gives (52.6 ±1.0, 13.3 ±0.8). 

3.4.2. Velocity Dispersion Profile 

To examine the velocity dispersion as a function of radius, 
we divide the face of Fornax into nine annuli containing ap- 
proximately equal numbers (19-21 per annulus) of member 
stars. From the stars in each annulus, we estimate the intrinsic 
radial velocity dispersion using the maximum likelihood tech- 
nique described above; here we modify the procedure, how- 
ever, so that the estimated mean true velocity, (m), of all bins is 
fixed at the value obtained from the global, unbinned sample. 
Since Fornax has a measured elUpticity of 0.3 ± 0.01 (IH95), 
we estimate the profiles using elliptical as well as circular an- 
nuli. The resulting radial velocity dispersion profile estimate, 
dp(R), for each of the three levels of Fornax membership re- 
jection is shown in Figure 5. The profiles using circular annuli 
display more scatter than those with elUptical annuU, although 
in general all are consistent with a flat profile to the linait of 
the sampled region. 

In addition to raising the overall dispersion, relaxing the 
membership criteria for the N=182 and N=186 samples em- 
phasizes an upturn in the dispersion at the outermost annulus. 
This feature persists when varying both the shape and num- 
ber of annuU, and is not Ukely an artifact of the HRF apparent 
rotation signal (Figure 4a). To demonstrate this last point we 
consider the outermost circular annulus from the N=182 ve- 
locity dispersion profile (Figure 5b). This annulus contains 21 
stars, spanning projected radii between 37 - 67 arcmin, and 
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Fig. 4. — Rotation signal of Fornax. The difference between the radial velocity of Fornax members on either side of a line passing through Fornax's center is plotted as a function 
of the position angle of that line, a) computed using the measured heliocentric rest frame radial velocities, uncorrected for perspective-induced rotation (see Section 3.3). b) computed 
from GRF radial velocities obtained using the Fornax proper motion measurement of Piatek et al. (2002). c) computed from GRF radial velocities obtained using the Fornax proper 
motion measurement of Dinescu et al. (2004). d) computed from GRF radial velocities obtained under the assumption that Fomax does not rotate. 



has, in units of km s"', ((m),(Tp) = (53.0 ± 3.4, 15.0 ± 2.4). 
The eight stars to the northeast of the HRF apparent rota- 
tion axis (Section 3.3) have ((m), ct^) = (55.3 ±5.0, 13.9± 3.6), 
while the thirteen stars to the southwest have {{u),ap) = 
(48.1 ±4.1,14.8±3.0). The nine stars from the latter group 
that are located nearest the southwest corner of Figure 2b 
have ((m),(T^) = (47.5 ±5.4, 16.1 ± 3.9). In each of these sub- 
annular regions the local velocity dispersion is equal (within 
statistical uncertainties) to the velocity dispersion measured 
for the whole annulus. Thus the adoption of a fixed mean ve- 
locity over the entire profile has not significantly inflated the 
calculated dispersion even where the effects of a velocity gra- 
dient would be strongest. Rather, the dispersion measured at 
large radius is dominated by localized velocity scatter. The 
possible rise, or at the very least, the lack of a falloff, in the 
outer dispersion may place Fomax in contrast with the Draco 
and Ursa Minor dSphs, for which velocity dispersion profiles 
have recently become available (Wilkinson et al. 2004; see 



section 5.2). 

Lin & Dong (2005) point out that a perspective-induced 
HRF radial velocity gradient may produce a discrepancy be- 
tween HRF and GRF radial velocity dispersion, particularly at 
large radii. Overplotted without errorbars in Figure 5 are the 
GRF radial velocity dispersion profiles measured after apply- 
ing Equations 5 and 6 to place the individual HRF velocities 
in the GRF, using either of the existing Fornax proper motion 
measurements. In all annuh the GRF velocity dispersion lies 
well within the Icr uncertainty region of the HRF dispersion. 
We conclude that for the present data set, the HRF radial ve- 
locity dispersion profile is a suitable surrogate for the GRF 
profile. 

3.4.3. Bias in the velocity dispersion estimate? 

By adopting Equation 7 to estimate velocity dispersion we 
implicitly assume that the stellar velocities everywhere fol- 
low a Gaussian distribution. This cannot strictiy be correct, as 
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Fig. 5 . — Radial velocity dispersion as a function of angular radius for three levels of Fornax membership discrimination (see section 3.2). Filled squares and errorbars correspond 
to HRF radial velocity dispersion. Open triangles and open circles indicate the GRF radial velocity dispersion calculated using, respectively, the Piatek et al. (2002) and Dinescu et al. 
(2004) values for Fornax's proper motion. The plots in the top row are constructed using circular annuli, while those in the bottom row use elliptical annuli with e = 1 — bja = 0.3, 
semimajor axis a = R and PA=41 " . a) Calculated using only the 176 stars with velocities surviving all four iterations of the bi-weight rejection algorithm, b) Calculated using the 182 
stars with velocities surviving the first two rejection iterations, c) calculated using the 186 stars with velocities surviving the first rejection algorithm. Bins contain approxunately equal 
munbers of stars. 



tidal fields will strip high-velocity stars and internal interac- 
tions will alter the velocity distribution with time. If the true 
stellar distribution function (DF) is non-Gaussian the velocity 
dispersion estimate, (Jp{R), may deviate systematically from 
the true velocity dispersion profile, properly calculated as the 
velocity moment of the stellar DF. 

To investigate the bias likely to be present in our estimate 
a-p(R) with respect to the profile calculated from a model 
DF, we again perform Monte Carlo simulations. For a given 
model we generate 1000 artificial data sets, each comprising 
stellar radial velocities for 186 stars occupying the same sky 
positions as the stars in our observed data set. The radial ve- 
locity assigned to each star is drawn randomly from the ap- 
propriate DF, integrated at each projected radius. For indi- 
vidual radial velocity uncertainties we adopt the same values 
as calculated for the observed sample in section 2.3. We bin 



each artificial data set, using circular annuli of the same radii 
and size as those in the observed N=186 data set. Within 
each annulus, we then calculate the velocity dispersion esti- 
mate, a-p, and the associated confidence interval for each of 
the 1000 subsamples. We sum the results to obtain a function 
that gives the probability of measuring a given ap within that 
annulus, in the case that Fornax adheres to the adopted model 
DF. We then compare to the velocity dispersion calculated di- 
rectly from the model DF. 

One example of a plausible non-Gaussian DF is given 
by King (1962, 1966; Binney & Tremaine 1987, hereafter, 
BT87). The King model has radius and velocity scale param- 
eters and Vj, as well as a third parameter, Wq, that specifies 
the value of the central potential in units of Vj. In our simu- 
lations we adopt King models with Wo = 3.26 and Wo = 10.0. 
These provide reasonable single-component fits to the surface 
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brightness and flat velocity dispersion profiles, respectively. 
The latter model is also close to an isothermal sphere and so 
should have httle distortion with respect to a Gaussian DF. 
In both cases we set = 13. T (IH95) and leave v, as a free 
parameter. 

Figure 6 shows the resulting probability functions for ap 
within each annulus, and identifies the model projected veloc- 
ity dispersion at the radius of the annulus. The projected ve- 
locity dispersion for the nearly Gaussian Wq = 10.0 model lies 
near the center of the simulated ap distribution in each annu- 
lus. For Wo = 3.26, the simulated ap values are sUghtly biased 
in favor of overestimating the model dispersion. The discrep- 
ancy between the Wo = 3.26 model velocity dispersion and the 
median simulated ap ranges from 0.04v5 at the outermost an- 
nulus to O.OSvj at the innermost. Since Vj roughly equals the 
central velocity dispersion, this discrepancy amounts to ~ 1 
km s~' for a dSph-Uke system. 

We conclude that, despite the formal distinction between 
the estimated velocity dispersion profile, ap(R), and the pro- 
jected velocity dispersion profile calculated from a model DF, 
the former provides an unbiased estimate of the latter in the 
case of a Gaussian DF. With respect to non-Gaussian DFs, 
(TpiR) may introduce a bias, and one should exercise caution 
when comparing ap{R) to such models. 

4. ANALYSIS 

4.1. Equilibrium Models 

The classical analysis of dSph velocity data sets has re- 
lied on application of equihbrium models falling under the 
purview of the core-fitting technique (Richstone & Tremaine 
1986). Chief among these is the single-component King 
model, which parameterizes the stellar DF under characteris- 
tic assumptions of dynamic equihbrium, spherical symmetry, 
velocity isotropy, and a mass profile that is directly propor- 
tional to the luminous profile ("mass follows light"). While 
they may provide a reasonable approximation of dSph cores, 
none of these assumptions are easily justified over the ex- 
tended regions sampled by modem data sets. The Milky Way 
dSphs have measured elhpticities ranging from 0.13 to 0.56 
(IH95). The degree to which ongoing tidal interactions with 
the Milky Way cause departures from equilibrium is con- 
troversial and poorly constrained. The presence of veloc- 
ity anisotropy results in a well-known degeneracy with mass. 
Setting aside modifications to Newtonian gravity, on no other 
galactic scales does mass follow light (BT87; Kormendy & 
Freeman 2004). 

The velocity dispersion profile we measure for Fornax pro- 
vides compelling evidence that the classical analysis is insuffi- 
cient. A single-component King model that assumes mass fol- 
lows Ught has, adopting the photometrically-determined For- 
nax structural parameters of IH95, Wo = 3.26 and = 13.7' 
(Figure 7a-c). The artificial data sets described in Section 
3.4.3 indicate that, for Wo = 3.26, the distribution of simulated 
ap values shifts to lower velocity dispersion toward larger an- 
gular radius (Figure 6). The median simulated ap in the out- 
ermost armulus drops to one-third the median simulated ap in 
the innermost annulus. We do not see evidence of this behav- 
ior in the actual data, for which the measured velocity disper- 
sion profile remains approximately flat at all radii. In none of 
the three Fornax membership samples is the outermost mea- 
sured ap less than the irmermost measured ap. We use the 
simulated ap probabiUty function in each annulus to calculate 
a negligible probability of measuring a flat velocity dispersion 
profile given the mass-foUows-light model. Figure 7d plots, 
for each annulus in the velocity dispersion profile, the prob- 



ability (from the simulated ap distribution) that ap is at least 
as large as the median simulated ap in the innermost annulus. 
The probabihty drops from 4% for the annulus at R=25.1' to 
2% at R=28.8', then falls by an order of magnitude at each of 
the two remaining annuli. We thus find a general failure of 
the single-component mass-foUows-light King model to re- 
produce the flat behavior of the observed velocity disperion 
profile. 

It is clear that at least one of the classical assumptions is 
invalid in Fornax. Perhaps the most readily discarded is mass 
follows Ught; indeed, a flat velocity dispersion profile may 
arise if the stars orbit inside a dark matter halo with core ra- 
dius larger than that of the visible component. We explore 
this scenario using two-component King models that continue 
to assume spherical symmetry, dynamic equilibrium, and ve- 
locity isotropy. These models contain the additional assump- 
tion of energy equipartition in the core region, which does not 
readily pertain to coUisionless systems such as dSphs. We 
therefore adopt an approach similar to that of Pryor & Kor- 
mendy ( 1 990), who used two-component King models merely 
as tools for exploring possible dark matter distributions in the 
Draco and Ursa Minor dSphs^. If £ = -v^/2-(/) is the total en- 
ergy per unit mass, and (p is the potential per unit mass, then 
for the isotropic case each component / has energy distribu- 
tion function (King 1966; Pryor & Kormendy 1990) 

fiiE) a g-wf^^-"'"'")/''' - 1 . (11) 

Two additional parameters join the familiar r,, v,, and W). 
For luminous (subscript l) and dark (subscript d) components. 
Pod/ Pol and jio/ P-l specify the ratios of central densities and 
dimensionless "masses," respectively. The /x, have a physi- 
cal interpretation when multi-component models are applied 
to stellar mass classes in globular clusters (e.g., Gunn & Grif- 
fin 1979; Da Costa & Freeman 1976). There, for mass class 
i, Hi = nii/rn, where rn is the sum of p'q poinii over all mass 
classes. For our purposes, jio / Pl determines the ratio of core 
radii, rco/rcL, given energy equipartition in the core. The 
"core radius," Tc is defined as the radius at which the pro- 
jected density of component ; falls to half its central value. 

We subject each model first to constraints set by the struc- 
tural parameters derived from the img/e-component King fit 
of IH95. We adopt r^L = 390 ± 36 pc (updated for a Fornax 
distance of 138 kpc) and So = 15.7±5.1 L0 pc~^ in the V 
band. The model surface brightness profile, Yj{R), is scaled 
by the product VsPoiM/L]'^, where po = Pod + Pol and [M/L]l 
is the V-band mass-to-hght ratio of the luminous component, 
assumed to be independent of radius. For a given model, we 
assign the value of r, to be that which places at the IH95 
value. We then assign /5o[M/L]^' the value that recovers the 
1H95 value for the central surface brightness. 

The models are next constrained by the available veloc- 
ity dispersion and surface brightness profiles. The values of 

and p()[M/L]J^ set the velocity scale according to 9vJ = 
AnGr^PQ (King 1966). For a given value of Wo and an adopted 
[M/L]l, we determine the {pod/ POL.fcD/ rcL) pair that pro- 
vides the best fit to the IH95 surface brightness profile while 

^ For an approach that does not assume energy coupUng between the lu- 
minous and dark matter, see Wilkinson et al. (2002), who use self-consistent 
stellar DFs to describe stars acting as tracers in a dark plus luminous potential. 
WiUdnson et al. and Pryor & Kormendy also allow for velocity anisotropy, 
whereas we consider only idealized, isotropic DFs. Radial velocity samples 
of the present size are only marginally able to address issues of anisotropy 
(Wilkinson et al. 2002), and we reserve a more comprehensive analysis for 
future work with larger data sets. 
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Fig. 6. — Test for bias in the velocity dispersion estimate, a^,, witli respect to the projected velocity dispersion calculated directly from a model distribution function. Lefl: 
Nine panels represent the annuli (with specified angular radius) of velocity dispersion profiles calculated from simulated data. Curves indicate the probability of measuring a velocity 
dispersion ap in each annulus, given 1000 simulated data sets drawn from a King model with Wq = 3.26 (solid) or Wq = 10.0 (dotted). The vertical line identifies the velocity dispersion 
calculated directly from the corresponding model DF. Right: Open squares indicate the median simulated d-p in each annulus. Errorbars enclose the 68% of simulated d-p values nearest 
the median. Curves represent the projected velocity dispersion calculated directly from the model DF. 
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Fig. 7. — (a)-(c) The projected velocity dispersion from a mass-follows-liglit King 
model (Wo = 3.26, = 13.7'; IH95) is drawn as a dashed curve over the measured ct^ (/J) 
(using circular annuli), for each of three Fornax membership samples. The solid ctffve in 
(b) is the velocity dispersion estimator of Wang et al. (2005; see section 4.2), which was 
calculated using this Fornax sample, (d) Plotted for each annulus is the probability, fmm 
simulated data drawn from the mass-follows-light King DF, of measuring a velocity 
dispersion at least as large as the velocity dispersion of the innermost annulus. 



having a central velocity dispersion equal to the global ve- 
locity dispersion of the Fornax sample. Of these models we 
consider those with velocity dispersion profiles remaining flat 
over the sampled Fornax region to provide the best overall 
agreement with the data. Thus we are approximating the ob- 
served velocity dispersion profiles of Figure 5 as perfectly flat 
and ignoring any bias in apiR) with respect to the model ve- 
locity dispersion profile. We find that the favored models tend 
to have large Wq values, for which such bias is expected to be 
minimal (Section 3.4.3). 



Models representing [M/L]i= 1,2,3 and a range of Wq are 
summarized in Table 6. The first column gives the number of 
stars considered to be members in the velocity sample. The 
next six columns list the adopted \M/L]i and model parame- 
ters. The eighth column gives per degree of freedom with 
respect to the velocity dispersion profile. We do not use the 
test to determine a "best-fit" model, but merely to indi- 
cate the degree to which the considered models are consistent 
with a flat dispersion profile. The final three columns list the 
derived quantities of interest: the central dark matter density, 
total mass, and overall V-band M/L. Projected velocity disper- 
sion and surface brightness profiles for these models are plot- 
ted in Figure 8. For simplicity we include only those models 
used with the N=182 Fornax sample. 

The available Fornax velocity data place several broad con- 
straints on the two-component models. First, poo must be 
of order poL or larger in order to recover the observed cen- 
tral velocity dispersion. This is best illustrated by the "two"- 
component models with poo / pql = (Figure 8a). These mod- 
els contain no dark matter, and are therefore equivalent to the 
single-component, Wq = 3.3 King model fit of IH95. Even if 
[M/L]l = 3, models lacking a dark component underpredict 
the central velocity dispersion by a factor of two, and fare 
much worse at larger radii. Thus it is difficult to explain the 
velocities of even the most central stars without invoking dark 
matter Recognizing that p^i cx [M/L]i, the models able to 
reproduce the data have central dark matter densities between 
0.04 - 0.10 M0 pc~^. This is similar to the model-independent 
lower Hmit of poo > 0.05 Mq pc"^, derived by Pryor & Kor- 
mendy for Draco and Ursa Minor 

Second, the models able to reproduce the data span a sur- 
prisingly narrow range in size of the dark halo, with 2 < 
rco/rcL < 3. The observed surface brightness profile helps 
set a Wo-dependent lower limit on rco/rcL- For models with 
Wo > 3.3 (i.e., models for which the central potential is deeper 
than a single-component fit to the luminous material would 
suggest), the shape of the luminous density profile is sensitive 
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Fig. 8. — Projected velocity dispersion (top row of panels) and surface brightness (bottom row) profiles calculated from two-component King models. From left to rigiit, tiie models 
have increasing values for the central potential parameter, Wo. Included in the velocity plots is the observed profile from the N-I82 Fornax sample. Each set of three curves corresponds 
to [M/L]l = 1, 2, 3; <Tp(R) falls off faster for larger [M/L]/,. The set of models at lower velocity dispersion in panel (a) have no dark matter, with Pod/ Pol = 0. The thick curve in the 
surface brightness plots is a single-component model having the IH95 structural parameters. The dotted curve in the lower panel of (d) represents a model with Wo = 9 wiircn/rcL = 1-1, 
and illustrates the effect of a small core radius ratio on the surface brightness profile. 



to the dark matter potential. If rco jfcL is sufficiently greater 
than unity, the dark matter density is constant over a few Vcl, 
and the luminous profile retains its shape even when increas- 
ing Wo or Pod- To demonstrate the disruptive effect of low 
rco/rcL on the surface brightness profile, the dashed line in 
the bottom panel of Figure 8d is the surface brightness profile 
for a model with r^D /''cL = 1 • 1 ■ 

The apparent upper limit on rco / may be a consequence 
of the assumption that the luminous and dark components are 
dynamically coupled. Under Equation 11, the rate at which 
Pil po, the density fraction of component i, decreases with ra- 
dius is determined by the value of /x,Wb. The density frac- 
tion for a component with large piWo declines sharply. For 
large r^o / rd, a model must have ^.lWq sufficiently larger than 
/xdWq. However, in the limit of small Pod/pol, Pl 1 as 
Pl/ Pd oo. Compared to the much less luminous Draco 
and Ursa Minor, Fornax favors models with smaller Pod/pol 



at smaller Wq. This tends to suppress plWq, thereby prevent- 
ing the luminous density profile from becoming much steeper 
than the dark matter profile. 

Finally, the flat velocity dispersion profile of Fornax favors 
models with large Wq. As Wq increases. Equation 1 1 tends 
toward the distribution function of a constant velocity disper- 
sion, isothermal sphere. The flattening of the resulting ve- 
locity dispersion profile is evident in Figure 8. Models with 
Wo > 7.0 remain sufficiently flat over the observed region, for 
any [M/L]l < 3, and so provide the best overall agreement 
with the data. Models having still larger Wq are not ruled out 
by the velocity data, although they eventually require larger 
values of [M/L]l. They also become unphysical, as M oc r for 
very large Wq. 

Mass and V-band [M/L] profiles for the most suitable Wq > 
7.0 and Wq = 9.0 models are shown in Figure 9. The total 
masses derived from these models fall in the range 4-18x10^ 
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Fig. 9. — Mass and V-band MjL profiles for two-component King models from 
Table 6 with Wq > 7.0. (a) gives the cumulative (luminous plus dark) mass profiles, (b) 
gives the cumulative mass-to-light ratio as a function of radius. The dashed curves are 
nonparametric estimates (Section 4.2), 



M©, or 3 —7 X 10 M© if we integrate the density profile only 
over the observed region r < 2 kpc. These are one to two or- 
ders of magnitude larger than previous mass estimates based 
on only central velocity dispersion data (see M91 and refer- 
ences therein). This dramatic increase is a consequence of 
the flat velocity dispersion profile. The two-component King 
models suggest that if this flatness arises from stars moving 
isotropically inside an extended dark halo, then Fornax is very 
dark even over the observed region R <\° , with M /L perhaps 
10 to 40 times larger than that of the luminous component. 

4.2. Non-Parametric Mass Estimation 

In Paper 1 we introduce a non-parametric method for esti- 
mating mass distributions from photometric and radial veloc- 
ity data. We assume spherical symmetry, velocity isotropy, 
and dynamic equilibrium. We do not assume mass follows 
light, nor do we adopt a parametric form for the stellar DF. 
Rather, we use the IH95 star count data and each velocity data 
point to estimate deprojected profiles for the stellar density, 
/(r), and squared velocity dispersion, /x(r). These relate to 
the underlying mass, M(r), via the Jeans equations (Equation 
4-55 of BT87). In addition to shedding some of the classi- 
cal assumptions, this technique offers the benefit of avoiding 
the problems inherent to binning a radial velocity data set (for 
a King analysis that avoids binning via maximum likelihood 
techniques, see Oh & Lin 1992). 

Briefly, M(r) is approximated as a spUne of the form 

m 

M{r) = Y.I3i{[r-rMP, (12) 

i=l 

in which the notation [x]+ indicates the greater of x or zero. 
The values (3i,...,f3m depend on /(r) and /i(r) and are es- 
timated using the available data and by imposing general 
shape restrictions on M(r) (e.g., M{r) is non-decreasing and 
M{r)r=o = 0). Here we add one further shape restriction to 
those described in detail in Paper I. Specifically, we require 
the mass density to be a decreasing function of radius (see 
Appendix for details). This gives a smoother M(r), eliminat- 



ing the plateau features present in the original Fornax estima- 
tion (see Figure 12 of the Appendix). Simulations indicate 
that a strong positive bias in the mass estimate arises beyond 
a radius enclosing ^95% of the measured stars. Inside this 
radius (^ 1.5 kpc for the present data set), the nonparametric 
technique gives M{r) accurate to within 20% when operating 
on data sets containing 1000 or more stellar velocities (see 
Figure 7 of Paper I). 

We apply the nonparametric technique to the present data 
set with the caveat that the uncertainty in M(r) will be at least 
a factor of two. The solid line in Figure 7b is th e non paramet- 
ric estimate of the velocity dispersion profile, ^//i(r), obtained 
using the N=182 Fornax sample and the star count data of 
IH95. The corresponding M(r) is given by the dashed line in 
the top panel of Figure 9. The nonparametric mass estimate 
has a larger central value and shallower slope than the two- 
component King models, but displays a similar mass and be- 
havior at radii larger than ~ 1 kpc. The bottom panel of Figure 
9 indicates thatM(r) rises less steeply than the luminosity pro- 
file until approximately the Fornax core radius, where the en- 
closed M /L reaches a minimum value of ^ 2 [M/L] © . Outside 
the core, dark matter dominates, reaching M/L ^ 15[M/L]© 
before the estimation terminates at a radius of 1.5 kpc. Larger 
and more extended data sets wiU be of great value in taking 
full advantage of this technique. 

4.3. External Tides 

Each of the mass models, as well as the nonparametric mass 
estimation technique we have apphed to the kinematics of 
Fornax, assumes that Fornax is in a state close to dynami- 
cal equilibrium. This may not be valid if the dSph steUar 
component is tidally heated as the galaxy orbits within the 
Milky Way potential. Claims of member stars projected at 
distances well beyond the nominal tidal radius (as determined 
by a single-component King model fit to photometric data) 
of some dSphs have been cited as possible evidence for tidal 
influence (IH95; Martinez-Delgado et al. 2001; Palma et al. 
2003; Munoz et al. 2005). Further, flie Sagittarius dSph, at a 
Milky Way distance of ~ 16 kpc, is clearly undergoing tidal 
disruption (Ibata et al. 1994; Mateo et al. 1996; Ibata et al. 
2001; Majewski et al. 2003), and so presents at least one case 
in which tides dominate. 

Various n-body simulations have addressed the degree to 
which the Milky Way's tidal influence on its satellites might 
alter their kinematics and derived masses. Oh et al. (1995) 
simulate the evolution of dSphs over several perigalacticon 
passages and conclude that even a tidaUy disrupted, unbound 
dSph steUar population may exhibit a velocity dispersion not 
significantly different from its pre-disruption value. Piatek 
and Pryor (1995) simulate single perigalacticon passages and 
add that even when a strong tidal encounter modifies the struc- 
ture and internal kinematics of a dSph, the core is least af- 
fected and the central M/L derived from the equilibrium as- 
sumption is virtually unchanged. The simulations of Klessen 
& Kroupa (1998) and Klessen & Zhao (2002) show that an 
unbound tidal renmant projected along the line of sight may 
display some of the kinematic and morphological features of 
dSphs; however, Klessen et al. (2003) later argue that the nar- 
row horizontal branch observed in Draco rules out the ubiq- 
uity of this scenario. 

The kinematic data set presented here gives an opportu- 
nity to examine certain predictions of tidal disruption mod- 
els. Along a disrupting satellite's orbit, stars nearest the par- 
ent system begin to lead the sateUite's center of mass as they 
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Fig. 10. — Orientation of the apparent rotation signal and published proper motion 
of Fornax. The solid line through (0,0) is the Fornax morphological major axis {PA = 
41° ± 1° (IH95)), and the dotted lines through the origin enclose the pi obable axes 
of the GRF apparent rotation signal (see section 3.3). The receding hemisphere is to 
the northeast. Arrows indicate the dnection and relative magnitude of Fornax's proper 
motion in the Milky Way rest frame, as measured independently by Ratek et al. (2002, 
solid arrow), and Dinescu et al. (2004, dashed arrow). The associated proper motion 
uncertainties are mapped conservatively, encompassing all directions allowed by the 
uncertainties quoted for the two components of proper motion. 



become unbound. Since the satellite's own gravity continues 
to pull on these stars in the direction opposite their motion, 
they lose energy in the reference frame of the parent system. 
The opposite holds true for the satellite's stars farthest from 
the parent system: as they become unbound, they trail the 
satellite's center of mass and gain energy in the parent's ref- 
erence frame, since the tug from the satellite is now in the 
same direction as their motion. The result is elongation along 
the satellite's orbit, and apparent rotation of the satellite about 
its minor axis as observed from the parent system (Piatek & 
Pryor 1995; Oh et al. 1995). Thus we might expect a disrupt- 
ing satellite to display two observables: a GRF rotation signal 
with axis of apparent rotation perpendicular to the morpho- 
logical major axis, and a proper motion parallel to the mor- 
phological major axis. 

For Fornax, the GRF rotation signal detected using the Pi- 
atek et al. proper motion (Section 3.3) has a rotation axis 
oriented nearly perpendicular to the galaxy's morphological 
major axis (Figure 10). The arrows in Figure 10 show the 
directions and relative magnitudes of the Piatek et al. and Di- 
nescu et al. proper motion measurements of Fornax. The Pi- 
atek et al. proper motion vector is clearly not aligned with the 
major axis, contrary to predictions from the models of tidal 
interaction. In contrast, the Dinescu et al. proper motion is 
more nearly parallel to the major axis. However, the ~ 1 .2 km 
s~^ GRF rotation implied by the Dinescu et al. proper motion 
is too slow to have statistical significance. Thus the two ob- 
servables indicative of tidal disruption are not simultaneously 
present in the existing kinematic data. 

If we ignore the orientation of rotation and proper motion 
and instead simply examine the major-axis velocity trend, we 
find no evidence for a tidally-induced velocity gradient along 
this axis. For a star with GRF radial velocity v, let D be the 
angular distance between the star and Fornax's minor axis 
(i.e., distance along the major axis), and let v,y, be the bulk 



Fig. 1 1 . — GRF mean radial velocity along Foi'nax's moi'phological major axis. 
GRF velocities are calculated using Fornax proper motion measurements by Piatek et al. 
(2002; panel a) and Dinescu et al. (2004; panel b). Solid lines represent the best-fitting 
radial velocity gradient, assuming cylindrical solid-body rotation. 



GRF radial velocity of Fornax. We model dv/dD, the GRF 
radial velocity gradient along the major axis, according to 
v{D) = Vsys + dv/dD. This assumes that any apparent rotation 
resembles that of a cylindrical solid body. We then use the 
unbinned GRF velocity data to solve for dv/dD via linear 
regression. Using the Piatek et al. proper motion, we find 
dv/dD ~ 0.1 km s~' arcmin"^ (3 km s~^ kpc~^), and the Di- 
nescu et al. proper motion gives a shallower dv/dD 0.03 
km s"^ arcmin"^ (0.8 km s"' kpc"'). These gradients are 
drawn over plots of the mean velocity along the major axis 
in Figure 1 1 . The unbinned data display significant scatter 
about either gradient, with each fit having x^/dof ~ 30. This 
contrasts with the predictions of tidal disruption models (see. 
Figure 7 of Piatek & Pryor 1995; Figure 5 of Klessen & Zhao, 
2002) in which tides produce ordered, monotonic, and typ- 
ically steeper gradients. The dispersion about either of the 
best-fit gradients is ~ 12 km s"\ identical to the overall sam- 
ple velocity dispersion. This indicates that the Fornax kine- 
matics are dominated by random motions, rather than the or- 
dered, streaming motions indicative of tidal disruption. 

It should be noted that tides are not the only mechanism by 
which a dSph might be altered from a state of dynamic equi- 
librium. During a search for extratidal structure, Coleman et 
al. (2004; 2005) discovered two lobes along the Fornax mi- 
nor axis and aligned with two shell-like features. They inter- 
pret these as signs of a recent merger event, rather than tidal 
tails. If a recent merger is confirmed, then there are local- 
ized regions within Fornax that have not had time to virialize, 
making irrelevant the concept of a virial tidal radius. To what 
extent such mergers may be pervasive in the dSph population 
remains highly uncertain, although Kleyna et al. (2003; 2004) 
have detected kinematically distinct substructure in Ursa Mi- 
nor and Sextans. In addition, Tolstoy et al. (2004) have 
found evidence for two populations of ancient stars with dif- 
fering metallicity and velocity distributions within the Sculp- 
tor dSph. It is clear that dSphs have more complicated his- 
tories than once thought, and many more stellar spectra are 
required to identify distinct components. 
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5. DISCUSSION 
5.1. Astrophysical Velocity Variability 

We have assumed thus far that the measured radial veloc- 
ities result exclusively from the underlying gravitational po- 
tential and kinematic status of Fornax. This is not necessarily 
true. Binary orbital motion may add a random component to 
any single-epoch velocity measurement, as may bulk stellar 
atmospheric motions. 

Including the M91 data, the Fornax data set now contains 
multi-epoch velocity measurements for 20 stars. Having no 
more than five, and in most cases only two, distinct measure- 
ments for any one of these stars, we cannot deduce binary 
parameters. Instead we identify binary candidates as those 
stars exhibiting velocity variability exceeding that which we 
would expect from the formal measurement errors. For each 
multiply-observed star, we calculate the Xofa obtained from 
the velocity measurements and their formal errors (see Tables 
3 and 4), as well as the probability, p(x^,f,^), that > xlhs 
given Gaussian random measurement errors. These probabil- 
ities should follow a uniform distribution between 0.0 and 1 .0 
if the stars are non-variable. For our sample of 20 multiply- 
observed stars, this would predict roughly two stars having 
pix^b.-) falling within each probability range 0.0-0.1, 0.1-0.2, 
0.9-1.0, but only 0.02 stars having p(xL) ^ 001, which 
is the probability threshold suggested by Olszewski et al. 
(1996, hereafter, OPA96) as indicative of binarity. While we 
do indeed find between 1-3 stars per each tenth in probability 
between 0.0 and 1.0, there are two stars — F-M20 and F2-9 — 
for which pixlj < 0.001. The excess over the expected num- 
ber of stars at very low pixlbs^ suggests that at least these two 
stars are exhibiting true velocity variability and are therefore 
binary candidates. This implies a binary "dicovery fraction" 
of 0.1 for the present Fornax sample. 

The actual Fornax binary frequency depends not only on the 
fraction of observed stars we can identify as binaries, but also 
on the efficiency with which we can identify binaries among 
our subset of multiply-observed stars. OPA96 perform sim- 
ulations over an expanse of binary orbital parameters in or- 
der to examine the discovery efficiency for a sample of 118 
stars with multiple velocity measurements in the Draco and 
Ursa Minor dSphs. We do not attempt to replicate their proce- 
dure regarding the present Fornax sample, primarily because 
OPA96 have superior statistics from multiple measurements. 
We wish to emphasize, however, the important conclusion 
from OPA96 that the presence of binary stars in a dSph ra- 
dial velocity sample ultimately has little effect on the derived 
velocity dispersion. The error in the velocity dispersion due 
to sampling uncertainty outweighs the error introduced by bi- 
naries. Simulations show that this result is a general feature 
of dSph-like velocity samples given measurement deviations 
similar to those in the Draco-UMi sample (OPA96 Table 7; 
see also Hargreaves et al. 1996, who reach a similar conclu- 
sion from binary simulations). 

A binary discovery fraction of 0. 1 suggests the presence of 
at least 20 unidentified binaries in our Fornax velocity sample 
of ~ 200 stars. We make a crude attempt to examine the ef- 
fect of unidentified binaries on the derived Fornax properties 
by observing the effects of removing the two known binaries 
from various subsamples of the velocity data set. First, if we 
calculate the velocity dispersion for the 20 stars with multi- 
ple velocity measurements, we obtain &muitipies = 15.8 ±2.7 
km s"'. If we remove the probable binaries F-M20 and F2- 
9 we calculate a slightly larger value, ^multiples = 16.4 ± 3.0 



km s~ . Considering the velocity dispersion profile, F-M20 
is in the innermost bin, and F2-9 is in the third bin. If we 
re-calculate the velocity dispersions in these bins after the re- 
moval of the candidate binary (both bins originally contain 20 
stars, so based on our binary discovery frequency, we expect 
at least 1 undetected binary to remain in each bin), we find 
that removal of F-M20 causes the dispersion estimate in the 
innermost bin to rise from 1 1 .3 ± 2.0 km s~' to 1 1 .5 ± 2. 1 km 
s~\ and the removal of F2-9 causes the dispersion estimate in 
the third bin to rise from 13.4 ± 2.4 km s"^ to 13.7 ± 2.6 km 
s~^. In both cases, inclusion of the probable binary has neg- 
ligible impact on the measured velocity dispersion. We can 
draw no strong conclusion from these tests, as the number of 
detected binaries is small and the discovery efficiency is un- 
known. Nevertheless, we find nothing to refute the conclusion 
of OPA96 that binary stars neghgibly inflate the measured ve- 
locity dispersion. 

A second possible source of radial velocity noise may come 
in the form of bulk motions in the atmospheres of the ob- 
served stars (Gunn & Griffin 1979). Pryor et al. (1988) find 
this velocity "jitter" to be as large as 4-8 km s"^ in globu- 
lar cluster red giants, although the effect appears to fall off 
rapidly in stars more than 0.5 mag dimmer than the tip of the 
giant branch. If we define the lip of the Fornax RGB to have 1 
~ 16.7 (Figure 1), then the region susceptible to velocity jit- 
ter includes 47 stars from the N=186 Fornax sample (all but 
one of these stars were originaUy selected for spectroscopic 
observation prior to the observation of the photometric data 
set presented in section 3.1). Of these 47, twelve have repeat 
velocity measurements that we may examine for variability. 
Only one of these, F-M20 — identified in the previous section 
as a binary candidate — has p(Xofo) ^ 0.01 . While the velocity 
variability of the dimmer binary candidate F2-9 is unlikely to 
be due to atmospheric jitter, the variability of F-M20 may be 
due in part to atmospheric motion. The distribution of p(Xofa) 
is otherwise uniform. The velocity dispersion calculated from 
the 47 brightest giants is 13.4 ± 1.5 km s~'. The velocity 
dispersion calculated from the remaining 139 stars from the 
N=186 sample is 12.9 ±0.8 km s"'. We conclude that atmo- 
spheric jitter does not have a significant impact on the Fornax 
radial velocities. 

5.2. Comparison with Draco and Ursa Minor 

Velocity dispersion profiles extending to the limits of the 
stellar distributions are now available for the dSphs Fornax, 
Ursa Minor, and Draco (Kleyna et al. 2002; Wilkinson et al. 
2004, hereafter, WK04;, Munoz et al. 2005). The behavior 
of the dispersion profiles at large radius is of great interest, 
capable not only of distinguishing between kinematic models, 
but also testing assumptions on which those models are based 
and searching for tidal influence on the kinematics. Although 
the profiles are generally flat, WK04 find sharp decreases in 
the velocity disperions of both Draco and Ursa Minor at the 
outermost point of both profiles. The observed drops are too 
sudden to be explained by isotropic King or Plummer profiles, 
although an abrupt change in the velocity anisotropy might 
explain the drop observed in Draco (see Mashchenko et al. 
2005), WK04 argue that anisotropy cannot by itself provide a 
plausible explanation for the more dramatic drop they witness 
in Ursa Minor. If the sharply declining dispersions in Draco 
and Ursa Minor are real features, they indicate an absence of 
tidal heating at the large radii at which they occur (Read et al. 
2005). Read et al. argue that dark matter halos of up to 10^ - 
lO"^ Mq are necessary to prevent tidal heating of the stars in 



16 



Walker et al. 



these dSphs. This would indicate similarity between Draco, 
Ursa Minor, and Fornax as described by the isotropic, two- 
component King models of section 4.1. It should be noted 
that Munoz et al. (2005) have recently reanalyzed the data 
of Wilkinson et al. and conclude that the presence of such a 
drop depends largely on the binning scheme and membership 
criteria employed. 

Regarding Fornax, we find no evidence for a decrease in 
the velocity dispersion at large radius for any birming scheme 
or reasonable membership criteria. Instead, we see some evi- 
dence for a mild increase at the outermost data point, partic- 
ularly as we allow stars from the wings of the observed ve- 
locity distribution into the Fornax membership. If this rise 
were attributable to tidal effects it would be puzzling that 
Fornax is susceptible to such external influence while Draco 
and Ursa Minor are not. Fornax has considerably larger lu- 
minous mass {Ly = 1.5 x IQ^Lq\ Mateo 1998) and lies at a 
greater distance from the Milky Way (the two proper mo- 
tion studies cited in this work estimate that the current dis- 
tance of ~ 140 kpc is near perigalacticon) than either Draco 
{Lv = 2.6 X 1O^L0, D = 82 ± 6 kpc) or Ursa Minor {Ly = 
2.9 X IQ^Lq, D = 66 ± 3 kpc). The degree to which the 
Milky Way influences the kinematics of its satellites remains 
an open and intriguing question. The remaining Galactic 
dSphs— Sculptor (Ly = 2.2 x IO^Lq, D = 79 ±4 kpc). Sextans 
(L = 5.0 X IO^Lq, Z) = 86 ± 4 kpc), Leo I (L = 4.8 x lO^L©, 
D = 250 ± 30 kpc), Leo II(L = 5.8x1 O^L© , D = 205 ± 1 2 kpc). 
Carina (L = 4.3 x IO^Lq, £) = 101 ± 5 kpc), and the recently 
discovered Ursa Major dSph (Ly ~ 4 x W^Lq, D ^ 100 kpc; 
Willman et al. 2005) — occupy a large region of parameter 
space. High-quahty velocity data sets are necessary to deter- 
mine any correlation between mass, orbital parameters, and 
the behavior of the dispersion profile at large radius. 

5.3. Summary and Conclusions 

We have presented new radial velocity measurements for 
156 (h-9 with uncertain membership status) stars belonging 
to the Fornax dSph. This increases the total number of For- 
nax stars with published velocities to 176 (-hIO). In order to 
test for rotation we have used existing Fornax proper motion 
measurements to place the heliocentric velocities in the For- 
nax rest frame. Adoption of the Piatek et al. (2002) proper 
motion results in a (marginally) statistically significant GRF 
rotation signal of ~ 2.5 km s~' about an axis at 1 12° ± 8°, near 
Fornax's minor axis. Adoption of the Dinescu et al. (2004) 
proper motion results in no statistically significant GRF rota- 
tion signal. Despite a favorable orientation with respect to the 
minor axis, the rotation signal stemming from the Piatek et al. 
proper motion is difficult to attribute to tidal influence, as the 
proper motion vector is perpendicular to Fornax's morpholog- 
ical major axis. Thus the two predictions from tidal disruption 
models — apparent rotation about the minor axis and elonga- 
tion along the satellite's orbit — are not simultaneously evident 
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in the present data. By examining localized velocity disper- 
sions and the velocity gradient along the major axis, we have 
demonstrated that the stellar kinematics of Fornax are domi- 
nated by random, rather than bulk rotational or streaming tidal 
motions. 

The Fornax radial velocity dispersion profile is generally 
flat. We have demonstrated the inability of single-component 
King models to account for the observed velocity disper- 
sion profile. We have applied isotropic, two-component King 
models consistent with the observed Fornax surface bright- 
ness profile, and found that models having similar central den- 
sities for dark and luminous matter are able to reproduce the 
flat observed velocity dispersion profile if the dark matter halo 
has a core of at least twice the size of the luminous material. 
Two-component models favored by the data have masses in 
the range M 10^ - 10^ M© . This would indicate a similarity 
between Fornax and the dSphs Draco and Ursa Minor, if ex- 
ternal tides indeed do not affect the stellar kinematics in these 
systems. In this case, dSphs are even more massive and dark- 
matter dominated than previously thought, which may help 
ameliorate the "accounting problem" faced by cold dark mat- 
ter models (Klypin et al. 1999; Moore et al. 1999; Stoehr et 
al. 2002). 

In any case, these results add to the mounting evidence that 
we must turn to more sophisticated analytical tools in order to 
explain the motions of dSph stars, particularly at large radii. 
We have discussed one such tool — a nonparametric statistical 
smoothing technique for estimating spherical masses directly 
from radial velocity data. An application of this method to the 
present Fornax sample yields a model-independent estimate 
of the mass profile. The result is consistent with the large 
Fornax masses suggested by the two-component models. This 
will become a powerful tool as dSph data sets continue to 
grow. We direct the interested reader to Paper I, in which this 
method is introduced in formal detail, and to the appendix to 
this work. 

We thank the staff at Las Campanas Observatory for sup- 
port during these observing runs. We are grateful to Slawomir 
Piatek and Tad Pryor for valuable comments and discussion. 
We also thank the referee, D.N.C. Lin, for comments and crit- 
icisms that improved this work. We thank Ewa Lokas for 
calling to our attention the fact that in the preprint version 
of this article, repeat observations of three stars were listed in 
Table 1 with unique stellar IDs. This has been corrected in 
the present version, and calculated quantities have been up- 
dated (results are not affected significantly). MM and MGW 
are supported in part by NSF grants AST-0507453, AST- 
0206081, and previously by NSF grant 0098661. EO is sup- 
ported in part by NSF grants AST-0098518, AST-0205790, 
and AST0507453. This work also received funding from the 
Horace H. Rackham School of Graduate Studies at the Uni- 
versity of Michigan, Ann Arbor (Rackham Interdisciplinary 
Collaboration Research Grant). 



Coleman, M.G., Da Costa, G.S., Bland-Hawthorn, J., Freeman, K.C. 2004, 
AJ, 127, 832 

Coleman, M.G., Da Costa, G.S., Bland-Hawthorn, J.B., Freeman, K.C. 2005, 

AJ, 129, 1443 
Da Costa, G.S., & Freeman, K.C. 1976, ApJ, 206, 128 
Dehnen, W., & Binney, J.J. 1998, MNRAS, 298, 387 
Dinescu, D.I., Keeney, B.A., Majewski, S.R., Girard, T.M. 2004, AJ, 128, 

687 

Evans, D.S. 1967, The Revision of the General Catalogue of Radial 
Velocities, lAUS, 30, 57 



Fornax Kinematics 



17 



Fleck, J.J., & Kuhn, J.R. 2003, ApJ, 592, 147 
Gunn, J.E., & Griffin, R.F. 1979, AJ, 84, 752 

Hargreaves, J.C., Gilmore, G., Annan, J.D. 1996, MNRAS, 279, 108 
Ibata, R.A., Gilmore, G., Irwin, MJ. 1995, MNRAS, 277, 781 
Ibata, R., Irwin, M., Lewis, G.F., Stolte, A. 2001, ApJ, 547L, 133 
Irwin, M., & Hatzidimitriou, D. 1995, MNRAS, 277, 1354 
King, I.R. 1962, AJ, 67, 471 
King, I.R. 1966, AJ, 71, 64 
Klessen, R.S. & Kroupa, P. 1998, ApJ, 498, 143 
Klessen, R.S. & Zhao, H. 2002, ApJ, 566, 838 
Klessen, R.S., Grebel, E.K., & Harbeck, D. 2003, ApJ, 589, 798 
Kleyna, J.T., WUMnson, M.I., Evans, N.W., Gilmore, G. 2001, ApJ, 563L, 
115 

Kleyna, J.T., Wilkinson, M.I., Evans, N.W., Gilmore, G., Frayn, C. 2002, 
MNRAS, 330, 792 

Kleyna, J.T., Wilkinson, M.I., Gilmore, G., Evans, N.W. 2003, ApJ, 588L, 21 
Kleyna, J.T., Wilkinson, M.I., Evans, N.W., Gilmore, G. 2004, MNRAS, 
354L, 66 

Kleyna, J.T., Wilkinson, M.I., Evans, N.W., Gilmore, G. 2005, ApJ, 630L, 
141 

Klypin, A., Kravtsov, A.V., Valenzuala, O., Prada, F 1999, ApJ, 522, 82 
Kormendy, J., & Freeman, K.C. 2004, Dark Matter in Galaxies, lAU Symp. 

220, eds. S.D. Ryder, D.J. Pisano, and K.C. Freeman, p. 377 
Kroupa, P 1997, New A, 2, 139 
Kuhn, J.R. 1993, ApJ, 409L, 13 
Lin, D.N.C. & Dong S. 2005, in progress 

Majewski, S.R., Skrutskie, M.F., Wemberg, M.D., Ostheimer, J.C.. 2003, 

ApJ, 599, 1082 

Martinez-Delgado, D., Alonso-Garcia, J., Aparicio, A. 2001, ApJ, 549, L63 
Mashchenko, S., Couchman, H.M.P, Sills, A. 2005, (astro-ph/0502I26) 
Mateo, M., Olszewski, E., Welch, D.L., Fischer, P., Kunkel, W. 1991, AJ, 
102, 914 

Mateo, M., Mirabel, N., Udalski, A., Szymanski, M., Kaluzny, J., Kubiak, 

M., Krzeminski, W., Stanek, K. 1996, ApJ, 458L, 13 
Mateo, M.L. 1998, ARA&A, 36, 435M 

Moore, B., Ghigna, S., Govemato, F, Lake, G., Quinn, T., Stadel, J., Tozzi, 

P. 1999, ApJ, 524L, 19M 
Munoz, R.R., Frinchaboy, P.M., Majewski, S.R., Kuhn, J.R., Chou, M.Y., 

Pahna, CP, Sohn, S.T., Patterson, R.J., Siegel, M.H. 2005, (astro- 

ph/0504035) 

Navarro, J.F, Frenk, C.S., White, S.D.M. 1997, ApJ, 490, 493 



Odenkirchen, M., Grebel, E.K., Harbeck, D., Dehnen, W., Rix, H-W., 
Newberg, H.J., Yanny, B., Holtzman, J., Brinkmann, J., Chen, B., and 13 
coathors, 2001, AJ, 122, 2538 

Oh, K.S., & Lin, D.N.C. 1992, JKAS, 25, 10 

Oh, K.S., Lin, D.N.C, Aarseth, S.J. 1995, ApJ, 442, 142 

Olszewski, E.W., Schommer, R.A., Suntzeff, N.B., Harris, H.C 1991, AJ, 
101, 5150 

Olszewski, E.W., Pryor, C, Armandroff, T.E. 1996, AJ, 111, 750 

Palma, C, Majewski, S.R., Siegel, M.H., Patterson, R.J., Ostheimer, J.C, 

Link, R. 2003, AJ 125, 1352 
Piatek, S. & Pryor, C. 1995, AJ, 109, 1071 

Piatek, S., Pryor, C, Olszewski, E.W., Harris, H.C, Mateo, M., Minniti, D., 

Monet, D.G., Morrison, H., Tinney, CG. 2002, AJ, 124, 3198 
Pryor, CP, Latham, D.W., Hazen, M.L. 1988, AJ, 96, 123 
Pryor, CP, & Kormendy, J. 1990, AJ, 100, 127 

Read, J.I., Wilkinson, M.I., Evans, N.W., Gilmore, G., Kleyna, J.T. 2005, 
(astro-ph/0505226) 

Rice, John A. 1995, Mathematical Statistics and Data Analysis, Duxbury 
Press 

Richstone, D.O., & Tremaine, S. 1986, AJ, 92, 72R 

Schechter, PL., Mateo, M., Saha, A. 1993, PASP 105, 1342 

Shectman, S.A. 1984, Proc. SPIE 445, Instrumentation in Astronomy V, e. A. 

Boksenberg & D.L. Crawford (Bellingham: SPIE), 128 
Stoehr, F, White, S.D.M., Tormen, G., Springel, V. 2002, MNRAS, 335L, 84 
Tolstoy, H., Irwin, M.J., Helmi, A., Battaglia, G., Jablonka, P., Hill, V., Venn, 
K.A., Shetrone, M.D., Letarte, B., Cole, A.A., Primas, F, Francois, P., 
Arimoto, N., Sadakane, K., Kaufer, A. 2004, ApJ, 617L, 119 
Tonry, J. & Davis, M. 1979, AJ, 84, 151 1 

Udry, S., Mayor, M., Maurice, E., Andersen, J., Imbert, M., Lindgren, H., 
Mermilliod, J.C, Bordstrom, B., Prevot, L. 1999, in Precise Stellar Radial 
Velocities, ed. J.B. Heamshaw & CD. Scarfe, ASP Conf. Ser., 185, lAU 
Coll., 170, 383 

Wang, X., Woodroofe, M., Walker, M.G., Mateo, M., Olszewski, E.W. (Paper 

I), 2005, ApJ, in press 
WiUdnson, M.I., Kleyna, J., Evans, N.W, & Gihnore, G. 2002, MNRAS, 330, 

778 

Wilkinson, M.I., Kleyna, J.T., Evans, N.W., Gihnore, G.F., Irwm, M.J., 

Grebel, E.K. 2004, ApJ, 61 IL, 21 
Willman, B., et al. 2005, ApJ, 626, 85 



APPENDIX 



Wang et. al. (2005) modelled M by a quadratic spline, 

m 



(1) 



!=1 



and reduced the estimation problem to a quadratic programming problem in which a quadratic function Qipi , • • • , (3m) is to be 
minimized subject to the constraints 



i=l 



and 



^A(o-'-m)>o 



(2) 



(3) 



i=l 



for j = I,- ■ ■ ,m-l. These constraints are necessary for a quadratic spline to be non-negative, non-decreasing, and bounded. 
Denoting the solution to the quadratic programming problem by /3i, • • • ,/3m, the resulting estimator, given by Equation 1, tracked 
the gross features of (the true) M quite well in simulations, but it often has flat stretches, implying regions of no estimated mass 
(Figure 12, solid Une). The latter feature can be eliminated (as in the dashed hne in Figure 12) by supposing that the mass density 
p{r) is non-decreasing in r. This assumption leads to the further constraints 



^A[2rM-'-d<0 



(4) 



for k = j-1 and k = j for j = I, - ■ ■ ,m. Estimation of M with the additional assumption that p is non-increasing then proceeds as 

in Paper I: the estimated M is given by (1), where now ,(}„ denote the solution to the quadratic programming problem with 
the constraints (4) in addition to (2) and (3) 
It is also possible to estimate the density, since 

1 
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where the prime denotes derivative. For M this becomes 

m , 



!=1 



Unfortunately, (5) becomes infinite as r ^ 0, a feature that is forced by the use of quadratic splines. This is not a fundamental 
problem as the total mass remains finite, nor is it unusual within the realm of widely used dynamical models (see p(r) for the 
NFW model (Navarro et al. 1997), or the classical isothermal sphere). 
To see how the additional constraints (4) arise, first observe that p is non-increasing if and only if 



For M of the form (1), 



1=1 

in the interval rj-i <r <rj for j = l,---,m. So, 

'•'^[^^W]=2X]A[2r,-i-r] 

!=1 

in the interval r^-i < r < for y = 1, • • • ,m. Viewed as a function on the interval [0,rm], the latter is a discontinuous, piecewise 
linear function. So, it wiU be non-positive everywhere if and only if it is non-positive at the endpoints of each interval rj-i <r<rj. 
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Fig. 12. — Nonparametric estimation of the Fornax mass. The solid curve is the original result published in Wang et al. 2005. The dashed curve is the estimate produced using the 
same data, under the additional constraint that the mass density is a non-increasing function of radius. 
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TABLE 1 

Heliocentric Radial velocity results for Fornax Target Stars 



Stnr 


1^2000 


''2000 


HJD 


T IT rintp 





P7t 


J 


— — 












(-2400000.0) 


of observation 


(arcmin) 


(degrees) 






(km s ) 




Fl-l 


U2:38;51.b 


-34:35:25.5 


48955.6 


29 Nov 1992 


14.5 


12U.2 


16.18 


2.21 


55.8 ±2.7 


Y 








52623.5 


15 Dec 2002 










64.5 ±2.0 




Fl-2 


02:38:52.3 


-34:44:55.3 


48955.7 


29 Nov 1992 


20.8 


143.7 


16.41 


1.64 


42.9 ±2.6 


Y 








48961.6 


05 Dec 1992 










45.2 ±2.2 




Fl-3 


02:38:59.9 


-34:45:26.9 


48955.7 


29 Nov 1992 


20.4 


148.1 


16.54 


1.72 


49.6 ±2.3 


Y 


Fl-4 


02:39:53.3 


-34:46:02.6 


48955.8 


29 Nov 1992 


17.9 


180.7 


16.41 


1.83 


44.7 ±2.6 


Y 








52621.6 


13 Dec 2002 










43.5 ±1.7 










52622.6 


14 Dec 2002 










43.5 ±1.7 




Fl-5 


02:38:03.0 


-34:16:39.1 


48956.6 
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25.3 
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16.74 


1.55 


45.5 ±2.7 


Y 
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9.5 
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Y 
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22.0 
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1.44 


48.5 ±2.8 


Y 
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02:39:28.1 
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5.4 
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1.63 


42.7 ±2.8 
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Fl-10 


02:39:52.3 
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02 Dec 1992 


9.7 
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1.48 


51.3±2.7 


Y 
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02:40:07.2 


-34:19:18.1 


48958.6 
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9.4 
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1.57 


55.1 ±2.6 


Y 
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11.1 


334.6 


16.72 


1.58 


41.5±2.8 
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48958.7 


02 Dec 1992 


30.7 
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16.78 
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Fl-16 


02:39:11.1 


-34:39:09.0 


48959.6 


03 Dec 1992 


13.9 
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16.74 


1.62 


60.6 ±2.7 
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02:39:36.8 


-34:45:29.7 


48959.7 
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17.6 


169.6 
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1.54 


67.7 ±2.6 
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02:39:40.0 


-34:43:00.9 


48959.7 


03 Dec 1992 


15.1 
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28.5 ±2.8 
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1.62 


35.5 ±2.8 
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1.47 
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11.2 
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1.43 
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12.9 
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1.41 


30.0 ±2.9 


Y 


Fl-28 


02:39:06.1 
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48960.7 


04 Dec 1992 


13.3 
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1.45 


45.2 ±2.7 


Y 
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48961.7 
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1.82 
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48961.8 
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48963.7 
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6.4 
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16.36 


1.79 
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Y 








49341.6 


20 Dec 1993 










43.8 ±2.0 










49649.6 


23 Oct 1994 










40.3 ±2.5 




F-M15 


02:39:54.0 


-34:34:24.3 


49333.6 


12 Dec 1993 


6.3 


183.1 


15.90 


1.96 


61.7±1.8 


Y 


F-M4 


02:40:01.8 


-34:27:48.1 


49333.6 


12 Dec 1993 


2.0 


280.1 


16.82 


1.63 


52.3 ±2.0 


Y 








49648.6 


22 Oct 1994 










53.4 ±2.6 




F18-1 


02:41:17.6 


-34:13:07.2 


49333.6 


12 Dec 1993 


23.1 


310.4 


16.98 


1.50 


50.0 ±2.0 


Y 


F15-1 


02:42:10.0 


-34:18:16.5 


49333.7 


12 Dec 1993 


30.1 


289.0 


17.09 


1.46 


59.4 ±2.3 


Y 


F15-2 


02:42:51.4 


-34:19:54.3 


49333.7 


12 Dec 1993 


37.8 


282.4 


16.78 


1.59 


60.9 ±2.3 


Y 


F18-2 


02:41:47.8 


-34:16:37.8 


49333.7 


12 Dec 1993 


26.5 


295.7 


16.86 


1.62 


56.7 ±2.3 


Y 


F18-3 


02:41:49.3 


-34:12:38.0 


49333.7 


12 Dec 1993 


28.7 


302.6 


17.07 


1.48 


66.4 ±2.6 


Y 


F-M20 


02:40:05.5 


-34:27:43.2 


49334.6 


13 Dec 1993 


2.8 


278.9 


16.01 


2.18 


63.7 ±2.3 


Y 


F-Ml 


02:39:39.6 


-34:19:52.2 


49334.6 


13 Dec 1993 


8.7 


17.5 


16.36 


1.95 


56.3 ±2.3 


Y 








49649.8 


23 Oct 1994 










51.8±2.8 


Y 


F19-2 


02:39:41.2 


-34:11:05.7 


49334.6 


13 Dec 1993 


17.2 


7.7 


16.87 


1.49 


56.1 ±2.2 


Y 


F19-3 


02:39:10.6 


-34:10:54.4 


49334.6 


13 Dec 1993 


19.3 


26.6 


16.77 


1.51 


20.1 ±2.3 


9" 


F20-1 


02:39:49.5 


-34:05:18.6 


49334.6 


13 Dec 1993 


22.8 


1.4 


16.97 


1.52 


57.2 ±2.7 


Y 


F17-1 


02:40:38.3 


-33:55:54.8 


49334.7 


13 Dec 1993 


33.6 


343.5 


16.78 


1.55 


134.0 ±2.7 


N" 


F20-2 


02:39:08.6 


-34:01:14.0 


49334.7 


13 Dec 1993 


28.4 


18.6 


16.91 


1.40 


59.6 ±2.4 


Y 


F20-3 


02:40:11.7 


-33:57:01.3 


49334.7 


13 Dec 1993 


31.4 


352.6 


17.02 


1.48 


44.2 ±2.6 


Y 


F20-4 


02:40:28.4 


-33:58:27.6 


49334.7 


13 Dec 1993 


30.6 


345.9 


16.99 


1.54 


64.6 ±2.5 


Y 


F17-2 


02:40:45.6 


-34:00:42.7 


49334.8 


13 Dec 1993 


29.6 


338.1 


16.76 


1.64 


49.5 ±2.7 


Y 


F-M2 


02:39:43.8 


-34:30:53.5 


49335.6 


14 Dec 1993 


3.2 


147.6 


16.26 


1.63 


71.8±2.6 


Y 


F13-1 


02:41:20.8 


-34:23:59.1 


49335.6 


14 Dec 1993 


18.7 


282.7 


16.80 


1.57 


66.2 ±2.1 


Y 


F13-2 


02:41:39.8 


-34:25:26.2 


49335.6 


14 Dec 1993 


22.3 


276.9 


19.43 


1.47 


47.2 ±2.3 


Y 


F13-3 


02:42:00.3 


-34:23:04.9 


49335.6 


14 Dec 1993 


26.9 


280.7 


16.78 


1.48 


68.2 ±2.7 


Y 


F13-4 


02:41:56.0 


-34:24:40.1 


49335.6 


14 Dec 1993 


25.7 


277.6 


16.87 


1.53 


61.7±2.2 


Y 


Fll-1 


02:41:28.6 


-34:44:19.0 


49335.7 


14 Dec 1993 


25.6 


230.7 


16.74 


1.67 


50.4 ±2.8 


Y 


Fll-2 


02:41:20.4 


-34:42:53.5 


49335.7 


14 Dec 1993 


23.4 


230.8 


16.88 


1.63 


46.0 ±2.5 


Y 


F13-5 


02:41:57.7 


-34:37:01.5 


49335.7 


14 Dec 1993 


27.3 


250.9 


16.85 


1.48 


46.3 ±2.2 


Y 


F14-1 


02:42:22.8 


-34:35:24.0 


49335.7 


14 Dec 1993 


31.8 


256.7 


16.90 


1.40 


60.7 ± 2.2 


Y 


Fll-3 


02:40:44.0 


-34:52:55.3 


49335.8 


14 Dec 1993 


27.0 


203.1 


16.82 


1.55 


53.9 ±2.2 


Y 


F8-1 


02:39:32.3 


-34:55:16.1 


49336.7 


15 Dec 1993 


27.4 


171.4 


16.98 


1.52 


54.2 ±2.6 


Y 


FlO-1 


02:41:21.4 


-35:00:01.8 


49337.6 


16 Dec 1993 


36.8 


209.8 


16.94 


1.58 


54.8 ±2.8 


Y 








52620.6 


12 Dec 2002 










53.0 ±2.2 




F9-2 


02:39:35.1 


-35:02:17.4 


49337.6 


16 Dec 1993 


34.3 


174.1 


17.04 


1.52 


45.2 ±2.2 


Y 


F9-3 


02:39:32.2 


-35:04:35.2 


49337.6 


16 Dec 1993 


36.7 


173.6 


16.94 


1.59 


52.1 ±2.3 


Y 


F9-4 


02:39:17.9 


-34:59:05.6 


49337.6 


16 Dec 1993 


31.7 


167.2 


16.85 


1.57 


55.2 ±2.0 


Y 


F9-5 


02:38:57.8 


-34:59:58.7 


49337.6 


16 Dec 1993 


33.7 


160.7 


16.77 


0.68 


65.3 ±2.4 




F6-1 


02:38:32.2 


-35:08:30.5 


49337.7 


16 Dec 1993 


43.6 


157.9 


16.92 


1.59 


32.9±2.1 


Y 


F6-3 


02:38:46.5 


-35:01:22.5 


49337.7 


16 Dec 1993 


35.9 


157.9 


16.83 


1.50 


86.2 ±2.8 


9" 


F6-4 


02:38:53.2 


-34:59:02.9 


49337.7 


16 Dec 1993 


33.2 


158.6 


17.03 


1.49 


47.8 ±2.4 


Y 


F9-6 


02:39:09.8 


-35:17:09.9 


49337.7 


16 Dec 1993 


49.8 


170.0 


17.14 


1.37 


68.3 ±2.5 


Y 








49651.9 


25 Oct 1994 










73.0 ±2.4 




F-M17 


02:39:46.5 


-34:25:52.7 


49338.5 


17 Dec 1993 


2.6 


27.8 


16.03 


1.53 


85.1 ±2.4 


?" 








49651.6 


25 Oct 1994 










86.3 ±2.4 










49652.5 


26 Oct 1994 










85.5 ±2.5 




F2-1 


02:37:48.4 


-34:36:21.9 


49338.6 


17 Dec 1993 


26.8 


108.0 


16.75 


1.58 


50.8 ±2.6 


Y 








49338.7 


17 Dec 1993 










50.1 ±2.3 




F22-1 


02:38:03.0 


-34:16:39.1 


49338.6 


17 Dec 1993 


25.3 


63.1 


16.74 


1.55 


50.5 ± 2.4 


Y 


F3-1 


02:36:32.5 


-34:28:13.8 


49338.6 


17 Dec 1993 


41.2 


90.3 


17.12 


1.47 


54.5 ±2.7 


Y 


F3-2 


02:36:53.9 


-34:25:52.6 


49338.6 
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07 Dec 1992 


68.5 ±0.9 












49333.5 


12 Dec 1993 


66.6 ±0.9 












49335.5 


14 Dec 1993 


68.1 ±1.0 












49337.5 


16 Dec 1993 


67.2 ±0.8 












49338.5 


17 Dec 1993 


68.0 ±0.8 












49339.5 


18 Dec 1993 


68.2 ±0.9 












49340.5 


19 Dec 1993 


66.7 ±0.8 












49341.5 


20 Dec 1993 


67.9±0.8 












49648.5 


22 Oct 1994 


66.1 ±0.9 












49649.5 


23 Oct 1994 


67.9±0.8 












49650.5 


24 Oct 1994 


68.2 ±0.8 












49651.5 


25 Oct 1994 


67.6 ±0.9 












49652.5 


26 Oct 1994 


68.0 ±0.9 
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TABLE 2 

Results for radial velocity standard stars 



CPD-432527 06:32:15.3 -43:31:14.3 +19.7 ±0.9' 



HD23214 03:42:09.1 -34:25:14.8 -A3±l.S^ 



HD43880 06:17:06.3 -34:44:13.1 443.6±2.42 



Twilight sky 



HD6655 01:05:18.0 -72:33:21.0 +19.5±0.3 



HD21581 03:28:54.8 
SA0217998 06:32:15.6 



-00:25:03.1 +154 ±1^ 
-43:31:13.4 +13. 1'* 



HD83516 09:38:02.9 -35:04:34.0 +43.5 ±0.2 



HD48381 06:41:43.0 -33:28:13.2 +40.5±0.2 



SAO201636 10:41:09.5 

HD2796 00:31:16.5 
HD103545 11:55:26.2 



-30:47:05.8 +262 ±1^ 

-16:47:43.6 -61 ± 1^ 
09:07:54.4 +180 ±1^ 



48954.9 


28 Nov 


1992 


19.7 ±0.7 


48955.9 


29 Nov 


1992 


20.4 ±0.7 


48957.8 


01 Dec 


1992 


19.8 ±0.8 


48958.9 


02 Dec 


1992 


20.1 ±0.8 


48959.9 


03 Dec 


1992 


20.2 ±0.8 


48960.9 


04 Dec 


1992 


18.1 ±0.8 


48961.9 


05 Dec 


1992 


20.0 ± 0.9 


48962.9 


06 Dec 


1992 


20.0 ±0.8 


48963.9 


07 Dec 


1992 


20.0 ±0.8 


49333.8 


12 Dec 


1993 


19.7 ±0.8 


49334.9 


13 Dec 


1993 


20.0 ± 0.8 


49335.9 


14 Dec 


1993 


22. 1 ± 0.9 


49335.9 


14 Dec 


1993 


20.0 ± 0.7 


49336.9 


15 Dec 


1993 


18.7±0.7 


49338.9 


17 Dec 


1993 


20.0 ±0.7 


49339.9 


18 Dec 


1993 


19.4 ±0.7 


49340 9 


19 Dec 


1993 


20 4 ± 7 


49341.9 


20 Dec 


1993 


20.4 ±0.8 


49648.9 


22 Oct 


1994 


19.7 ±0.7 


49649.9 


23 Oct 


1994 


19.7 ±0.7 


49650.9 


24 Oct 


1994 


20.6 ±0.7 


49651.9 


25 Oct 


1994 


19.1 ±0.8 


49652.9 


26 Oct 


1994 


18.7 ±0.9 


49653.9 


27 Oct 


1994 


18.9 ±0.8 


49333.6 


12 Dec 


1993 


-4.9 ±0.7 


49334 5 


13 Dec 


1993 


_5 3 -|_ Q 7 


49336^8 


15 Dec 


1993 


-5.3 ±0.1 


49648.9 


22 Oct 


1994 


-5.0 ±0.8 


49333.8 


12 Dec 


1993 


45.3 ±0.8 


49334.8 


13 Dec 


1993 


47.9 ±0.8 


49335.8 


14 Dec 


1993 


47.4 ±0.8 


49336.8 


15 Dec 


1993 


45.1 ±0.7 


49337.9 


16 Dec 


1993 


45.9 ±0.7 


jo.o 


i / uec 






49339.9 


18 Dec 


1993 


46.0 ±0.7 


49340.9 


19 Dec 


1993 


45.9 ±0.7 


49341.9 


20 Dec 


1993 


45.8 ±0.7 


49333 9 


12 Dec 


1993 


-1 4± 1 4 


49334.5 


13 Dec 


1993 


-0.6 ±1.0 


49335.5 


14 Dec 


1993 


0.7± 1.1 


49335 5 


14 Dec 


1993 


2 ± 1 


49337.5 


16 Dec 


1993 


-6.8 ±1.0 


49338.5 


17 Dec 


1993 


-0.3 ±1.0 


49339.5 


18 Dec 


1993 


-0.3 ±1.0 


49341.5 


20 Dec 


1993 


0.0 ±1.0 


52620.5 


12 Dec 


2002 


18.8±0.7 


52620.5 


12 Dec 


2002 


18.9±0.7 


52621.5 


13 Dec 


2002 


20.0 ±0.7 


52622.5 


14 Dec 


2002 


19.4 ±0.7 


52620.6 


12 Dec 


2002 


152.1 ±0.8 


52621.6 


13 Dec 


2002 


150.6 ±0.9 


52620.6 


12 Dec 


2002 


18.9±0.5 


52621.6 


13 Dec 


2002 


18.4 ±0.7 


52621.8 


13 Dec 


2002 


19.6 ±0.5 


52622 6 


14 Dec 


2002 


19 1 ±0 5 


52622^8 


14 Dec 


2002 


17.8±0.6 


52623.5 


15 Dec 


2002 


20.1 ±0.6 


52623 6 


15 Dec 


2002 


10 o±0 6 


52623.6 


15 Dec 


2002 


20.2 lb 0.6 


JZ.yjZ.J . o 


15 Dec 


2002 






iz uec 


9009 


Al, 1 -\-C\ A 


52621.8 


13 Dec 


2002 


42.4 ±0.4 


52622.9 


14 Dec 


2002 


42.9 ±0.5 


52622.9 


14 Dec 


2002 


42.6 ±0.4 


52622.9 


14 Dec 


2002 


43.9 ±0.6 


52623.8 


15 Dec 


2002 


40.2 ±0.6 


52621.8 


13 Dec 


2002 


40.2 ±0.5 


52623.8 


15 Dec 


2002 


40.3 ±0.6 


52621.8 


13 Dec 


2002 


264.3 ±0.9 


52623.5 


15 Dec 


2002 


-64.5 ±1.0 


52623.9 


15 Dec 


2002 


177.0 ±1.0 



+19.8 ±0.9 



-5.1 ±0.9 



+46.3 ±0.9 



-0.2 ±1.2 



+19.2 ±1.1 

+151.3 ±1.3 
+19.0±1.6 



+42.6 ±0.7 



+40.2 ±0.9 

+264.3 ±1.4 

-54.5 ±1.5 
+177.0±1.5 



^ published radial velocity 

* weighted mean measured radial velocity from the results presented in this table 
' reference: CORAVEL stody (Udry et al. 1999) 

^ reference: Olszewski et al. (1991) 
^ reference: Beers et al. (2000) 

* reference: Evans (1967) 
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TABLE 3 
Repeat velocity measurements. 



Star 




A. 




N 




(km s"^) 








(JPU-432b27 


19.8 ±U.8 


21.8 


0.532/ 


24 


Twffight Sky 


-0.2 ±1.0 


1.3 


0.9722 


7 


Fl-4 


43.7 ±1.8 


0.2 


0.9204 


3 


Fl-2 


44.2 ±2.4 


0.4 


0.5105 


2 


Fl-1 


61.4±2.3 


6.6 


0.0103 


2 


CPD-35°919 


4.1 ±0.7 


47.6 


0.0288 


32 


F-M17 


85.6 ±2.4 


0.1 


0.9364 


3 


F-M26 


42.0 ±2.2 


3.6 


0.3080 


3 


F-M4 


52.7 ±2.2 


0.1 


0.7369 


2 


F-M6 


34.0 ±2.2 


3.1 


0.0785 


2 


F-M8 


53.8 ±2.5 


2.6 


0.1101 


2 


FlO-1 


53.7 ±2.4 


0.3 


0.6086 


2 


F15-4 


54.8 ±1.9 


0.4 


0.5157 


2 


F17-4 


49.5 ±2.6 


7.3 


0.0069 


2 


F19-5 


46.9 ±2.8 


4.7 


0.0295 


2 


F2-1 


50.4 ±2.4 


0.0 


0.8516 


2 


F2-9 


50.9±2.1 


21.4 


0.0000 


2 


F24-2 


17.7 ±2.5 


2.5 


0.1137 


2 


F3-5 


55.4 ±2.6 


0.1 


0.7647 


2 


F31-3 


5.9 ±2.6 


3.2 


0.2063 


3 


F9-6 


70.7 ± 2.5 


1.9 


0.1734 


2 


HD 196983 


-9.3 ±0.7 


12.7 


0.6242 


16 


HD21581 


151.3 ±0.8 


1.5 


0.2151 


2 


HD2 19509 


67.8 ±0.9 


22.2 


0.3862 


22 


HD23214 


-5.1 ±0.7 


0.3 


0.9608 


4 


HD43880 


46.3 ±0.7 


14.6 


0.0685 


9 


HD48381 


40.2 ±0.5 


0.0 


0.8648 


2 


HD6655 


19.2 ±0.7 


1.9 


0.5916 


4 


HD83516 


42.6 ±0.5 


23.6 


0.0003 


6 


SA0217998 


19.0 ±1.0 


17.0 


0.0306 


9 



Probability that would have at least the measured value, given N independent measurements having the estimated uncertainties. 

TABLE 4 

Comparison with previously measured velocities. Entries with UT date prior to 
1992 were published in Mateo et al. (1991). All other entries are from the present 

study. 



Star 


UT Date 
of observation 


V 

km s-i 


km s-i 




P(x') 


F-Ml 


30 Nov 1989 
13 Dec 1993 
23 Oct 1994 


56. 7 ±2.9 
56.3 ±2.3 
51.8±2.8 


55.1 ±2.6 


2.0 


0.3/34 


F-MIO 


03 Dec 1989 

17 Nov 1990 

18 Dec 1993 


65.1 ±1.3 
63.8±1.8 
56.5 ±2.2 


63.1 ±1.6 


5.8 


0.0562 


F-M13 


03 Dec 1989 
25 Oct 1994 


92.1 ±2.0 

89.2 ±2.3 


90.9 ±2.1 


0.9 


0.3414 


F-M15 


04 Dec 1989 
12 Dec 1993 


57.9 ±2.0 
61.7±1.8 


60.0 ± 1.9 


2.0 


0.1579 


F-M17 


05 Dec 1989 
14 Nov 1990 
17 Dec 1993 

25 Oct 1994 

26 Oct 1994 


87.8±1.9 
86.7 ±2.2 
85.1 ±2.4 
86.3 ±2.4 
85.5 ±2.5 


86.5 ±2.2 


0.2 


0.9931 


F-M18 


05 Dec 1989 
05 Dec 1992 


54.3 ±2.4 
55.0 ±2.8 


54.8 ±2.6 


0.1 


0.8231 


F-M2 


01 Dec 1989 
14 Dec 1993 


71.6±1.5 
71.8 ±2.6 


71.6±1.8 


0.0 


0.9469 


F-M20 


15 Nov 1990 
13 Dec 1993 


38.9 ±2.0 
63.7 ±2.3 


49.6 ±2.1 


66.2 


0.0000 


F-M26 


18 Nov 1990 
07 Dec 1992 
20 Dec 1993 
23 Oct 1994 


40.3 ±1.6 
41.2±2.3 
43.8 ±2.0 
40.3 ±2.5 


41.4±1.9 


2.1 


0.54996 


F-M3 


01 Dec 1990 
24 Oct 1994 


59.1 ±2.8 
57.8 ±1.8 


58.2±2.1 


0.2 


0.6961 


F-M4 


01 Dec 1989 
15 Nov 1990 
12 Dec 1993 
22 Oct 1994 


60.0±3.6 
52.9±2.1 

52.3 ±2.0 

53.4 ±2.6 


53.6 ±2.3 


1.2 


0.7451 


F-M6 


02 Dec 1989 
16 Nov 1990 
18 Dec 1993 
26 Oct 1994 


34.9 ±1.9 
35.5 ±2.9 
31.9±2.0 
37.7 ±2.7 


34.5 ±2.2 


1.1 


0.7805 


F-M7 


02 Dec 1989 
24 Oct 1994 


42.3 ±2.5 
41.5 ±3.0 


42.0 ±2.7 


0.0 


0.8377 
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TABLE 4 



Comparison with previously measured velocities. Entries with UT date prior to 

1992 WERE published IN MATEO ET AL. (1991). ALL OTHER ENTRIES ARE FROM THE PRESENT 







study. 




F-M8 


02 Dec 1989 
04 Dec 1989 
19 Dee 1993 
27 Oct 1994 


54.4 ±2.5 54.6 ± 2.5 1.2 0.7644 
56.9 ±2.7 
51.4±2.3 
57.0 ±2.7 




TABLE 5 

ROTATION SIGNALS IN THE GALACTIC REST FRAME. 


Source 


Proper Motion ukP Rotation 
/*/ /ij, Speed Axis 
(mas century"' ) (mas century"' ) (kms"') (°) 


Signmcance 


Piatek et al. (2UU2) 
Dinescu et al. (2004) 


32±13 
-13 ±16 
-52 


33 ±13 2.5 ±0.4 112 ±8 
34±16 1.2±0.4 116±12 
+41 0.2 


10% 
74% 
100% 



TABLE 6 

Parameters and derived values from example two-component King models. 





LM/LJ/ 


Wo 


Pod/ POL 


rco/rcL 




V, 


X^/dof 


POD 




M/L" 




([M/L]q) 








(pc) 


(kms-') 




(Mq pc-3) 


(1O*M0) 


([M/L]q) 


FTF. 

i /D 




— TT — 

J.J 


rm 




— — 


J.J 


TTJ 


rm 


u.uu, u.uu 


1 


1 


176 


J 


3.3 


3.5 


1.7 


772 


17.5 


2.2 


0.067 


2.3; 2.2 


17 


16 


176 




5.0 


3.2 


2.1 


775 


17.0 


1.6 


0.061 


4.4; 3.7 


31 


27 


176 




7.0 


3.1 


2.2 


732 


15.8 


1.5 


0.059 


6.5; 4.1 


46 


30 


176 




9.0 


3.2 


2.3 


722 


15.7 


1.6 


0.060 


11.0; 4.4 


74 


31 


176 


2 


3.3 


0.0 




512 


7.8 


13 


0.31 


0.00; 0.00 


2 


2 


176 


2 


3.3 


1.4 


1.8 


729 


17.1 


2.5 


0.053 


2.0; 2.0 


15 


15 


176 


2 


5.0 


1.2 


2.3 


717 


16.1 


1.7 


0.046 


3.8; 3.4 


28 


26 


176 


2 


7.0 


1.2 


2.6 


683 


15.3 


1.6 


0.046 


5.8; 4.4 


42 


32 


176 


2 


9.0 


1.2 


2.5 


639 


14.4 


1.5 


0.046 


7.4; 4.3 


52 


32 


176 


3 


3.3 


0.0 




512 


9.6 


10 


0.47 


0.00; 0.00 


3 


3 


176 


3 


3.3 


0.8 


1.5 


632 


15.8 


3.4 


0.046 


1.2; 1.2 


12 


12 


176 


3 


5.0 


0.7 


2.0 


610 


14.8 


2.0 


0.040 


2.4; 2.3 


18 


18 


176 


3 


7.0 


0.6 


2.5 


589 


13.9 


1.7 


0.035 


4.1; 3.2 


27 


24 


176 


3 


9.0 


0.6 


3.1 


601 


14.1 


1.7 


0.035 


7.1; 4.8 


51 


36 


182 


1 


3.3 


4.3 


1.8 


852 


20.9 


2.3 


0.082 


3.8; 3.6 


27 


26 


182 


1 


5.0 


4.2 


1.9 


760 


18.5 


1.8 


0.080 


5.1; 4.2 


35 


30 


182 


1 


7.0 


4.1 


2.2 


766 


18.5 


1.7 


0.078 


9.3: 5.6 


67 


41 


182 


1 


9.0 


4.3 


2.1 


709 


17.4 


1.7 


0.082 


14.^; 5.9 


98 


35 


182 


2 


3.3 


1.8 


1.9 


791 


20.0 


2.5 


0.069 


3.2; 3.1 


24 


24 


182 


2 


5.0 


1.7 


2.1 


726 


18.1 


2.0 


0.065 


4.7; 4.2 


34 


30 


182 


2 


7.0 


1.6 


2.4 


699 


17.1 


1.7 


0.061 


7.2; 5.2 


52 


38 


182 


2 


9.0 


1.6 


2.6 


681 


16.6 


1.6 


0.061 


10.5; 5.8 


76 


42 


182 


3 


3.3 


1.1 


1.7 


681 


18.3 


3.3 


0.064 


2.0; 2.0 


15 


15 


182 


3 


5.0 


0.9 


2.1 


654 


16.8 


2.1 


0.052 


3.6; 3.3 


27 


24 


182 


3 


7.0 


0.8 


2.7 


644 


16.0 


1.7 


0.046 


6.1; 4.8 


45 


36 


182 


3 


9.0 


0.8 


2.9 


620 


15.5 


1.6 


0.046 


8.5; 5.4 


60 


39 


186 


1 


3.3 


5.1 


1.8 


857 


22.6 


2.1 


0.097 


4.4; 4.2 


32 


30 


186 


1 


5.0 


4.9 


2.0 


778 


20.2 


1.7 


0.093 


6.2; 5.0 


43 


36 


186 


1 


7.0 


4.8 


2.1 


750 


19.3 


1.5 


0.092 


10.0; 5.8 


70 


42 


186 


1 


9.0 


5.0 


2.1 


722 


18.9 


1.6 


0.095 


18.0; 5.9 


120; 41 


186 


2 


3.3 


2.1 


2.0 


835 


22.2 


2.2 


0.080 


4.2; 4.1 


32; 30 


186 


2 


5.0 


2.1 


2.1 


736 


19.6 


1.8 


0.080 


5.6; 4.8 


40; 36 


186 


2 


7.0 


2.0 


2.2 


694 


18.2 


1.6 


0.076 


8.0; 5.6 


56; 40 


186 


2 


9.0 


2.0 


2.4 


684 


17.9 


1.5 


0.076 


12.d; 6.2 


86; 44 


186 


3 


3.3 


1.3 


1.8 


725 


20.4 


2.7 


0.075 


2.8; 2.8 


21; 21 


186 


3 


5.0 


1.1 


2.3 


703 


18.9 


1.9 


0.063 


5.1; 4.6 


36; 36 


186 


3 


7.0 


1.1 


2.7 


687 


18.4 


1.7 


0.063 


8.5; 6.6 


63; 48 


186 


3 


9.0 


1.1 


2.7 


645 


17.3 


1.6 


0.063 


10.9; 6.6 


78; 48 



Number of member stars in Fornax sample 

Adopted mass-to-light ratio of the luminous component 

Mass of the dark component. The first value is the total mass of the dark component; the second value is the dark mass inside r< 2500 pc. 
V-band Mass-to-light ratio in solar units. The first value is the global M/L; the second value is M/L inside r< 2500 pc. 



